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Abstract 



We investigate the Ki Knn strong interaction decays. Using the ^Pq quark-pair- 
creation model to derive the basic parametrization, we discuss in detail how to obtain the 
various partial wave amplitudes into the possible quasi-two-body decay channels as well 
as their relative phases from the currently available experimental data. We obtain the 
Ki mixing angle to be Ok^ — 60°, in agreement with previous works. Our study can be 
applied to extract the information needed for the photon polarization determination of 
the radiative B ^ Kij decay. 
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1 Introduction 



1 Motivation for revisiting Ki-meson strong decays 

It has been proposed a method to measure the polarization of the photon in weak radia- 
tive decays of the S-meson by exploiting the decay B — )■ KTCTC'-f, with the system Kmr 
resonating into a Ki state [1,2]. We have recently extended this work [3] and shown 
that exploiting the full Dalitz plot for the Ktttt system could increase the sensitivity to 
the polarization determination inspired by the DDLR method [4]. For this purpose, it is 
important to have a good understanding of the strong i^i-decays. Indeed, it turned out 
that the B — )■ i^i (1270)7 channel, not considered in the original works [1,2], dominates 
over i^'i(1400) [5] while the pattern of the partial wave is especially complex for i^'i(1270). 

In the present paper, we give a detailed account and full discussion of the Ki hadronic 
decays. The amplitude of the process Ki — )■ Ktctt can be described by the basic quantity 
J: 

M{Ki^K7VTT)oce-J (1) 

where e is the polarization vector of the Ki in the Ki rest frame. The general framework 
for calculating J' is the quasi-two-body approximation: the process Ki — )■ Ktttt is decom- 
posed into two steps: 1) the decay of Ki — )■ vector isobar (V) + pseudoscalar (P); 2) the 
decay of the vector isobar {K* or p) into 2 pseudoscalars. Then, J7 is a sum of terms which 
are products of couplings and one isobar denominator ^. The explicit expressions for J' 
have been given in [3]. The decay properties of the intermediate isobars are well known. 
Here, we are then interested in evaluating the couplings describing the first step of the 
decay, Ki — )■ VP, and the relative signs or phases between the various channels. Concern- 
ing our motivation, it has appeared that the determination of the polarization parameter, 
called A^, of i? — )■ Kij depends essentially on the expression Im[n- (JT" x J'*)]/\J'\'^ [1^3]. 
This expression vanishes unless complex phases are present in J'. These are mainly pro- 
vided by the Breit-Wigner (BW) denominators of intermediate resonances (the so-called 
"isobars"), and possibly by complex phases of the couplings. It is found that such a 
quantity is very sensitive to the relative signs of the various channels in the strong decay, 
whence it is important to determine the signs or, possibly, complex phases of the cou- 
plings, which anyway can be observed in the various experiments provided one measures a 
sufficient number of angular distributions. One requires also specifically good knowledge 
of 1) D-waves; 2) off-shell extrapolation. 

2 Status of experimental study of Ki-mesons 

In principle, all the necessary hadronic parameters (i.e. Ki masses and partial decay 
widths, form factors and relative phases) can be determined from fits to the experimental 

^In the full expression for the weak process, e ■ J has still to be multiplied by a production amplitude 
and the corresponding Breit-Wigner denominator for the K\ . 
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data. However, at present moment we are far from being able to perform this with 
good accuracy, the experiments suffer from many drawbacks. The main reason is that 
the case under study cumulates many difficulties and complications, which have been 
underestimated in theoretical discussions: 

1. several possible partial waves {S,D) for the same channel 

2. three-body decay through multiple interfering channels 

3. broad parent Ki resonances 

4. effect of the large widths of isobars K* and p 

5. effect of a threshold (Kp) close to the i^'i(1270) resonance 

6. overlapping and mixing of two close states, fri(1270) and fri(1400) 

Perhaps not surprisingly for this particularly complicated case, one observes a rather 
confusing situation in experiments (e.g. contrary statements on total widths, on the 
KTT decay channel, etc.). This has been causing misunderstandings on these important 
observables. Then, one cannot simply use for example, the PDG entries, as is done usually. 
One should return to the original papers to understand what is actually measured, and 
this is not always easy. And also, we noticed some definite weak points. Up to now the 
most complete and accurate experimental analysis is the one by Daum et ai. It relies 
unhappily on the problematic phase space treatment of Nauenberg and Pais for decays 
to isobars. One notes also the absence of non resonant background in the i^-matrix, and 
the presence of unexplained "offset" phases. All this is explained in details in section 4.2. 
Finally, there is also lack of important informations like conventions of coupling signs 
or incomplete report of the parameters of the fit. Other experiments, which have been 
mentioned above, give precious complementary information but they are not able to solve 
all the problems, all the more since they are less accurate, and include less physical features 
in their fits (e.g. neglect of Z)-waves). Then, a sizable part of our work has consisted in 
an extensive discussion of the experimental analyses (mainly the one of Daum et al. and 
the one of Belle). 

3 The theoretical treatments of Ki-mesons 

The theoretical model has a first aim to give a physical understanding of the observed 
decay properties, which are far from trivial. In addition, it may serve to complement 
the experimental knowledge, where there are some lacks or weaknesses, and to help for 
future experimental analyses. Of course, there is no fundamental theoretical treatment 
of such processes. We have only phenomenological approaches at our disposal, mainly the 
one provided by the quark models. Approximate as it is, the quark model can be very 
precious to check the consistency of the present data and to orient the future studies 
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of i^i-decays. However, we have to keep in mind that it suffers from inherent, sizable 
and unknown uncertainties, which could limit the accuracy of the determination as 
mentioned earlier. To our knowledge, the best phenomenological model in the present 
case of ii'i-decays is the family of quark models, which are able to master the large set of 
hadronic states and their decays with a limited number of adjusted parameters ^. 

Admittedly, quark models are many, but one must distinguish between the potential 
models, and the decay models. The diversity is especially the one of potential mod- 
els, which intend to describe the spectroscopy of states. As concerns decays, there are 
not so many basic models. In fact, there are elementary emission models and quark- 
pair-creation models, both concerning two-body decays. This is why the quasi-two-body 
decay assumption is a natural step in the theoretical treatment of the three-body decay. 
The quark-pair-creation models have the advantage of unifying the whole of two-body 
and quasi-two-body decays. Among them, the ^Pq model (see section 3.2 and references 
therein) is particularly favoured as being the easiest to handle, and then the more ex- 
tensively tested, with a striking overall success over hundreds of decays. We use the 
^Po niodel with the important additional input of a damping factor, to account for off- 
shellness. Of course, as said before, quark models are inherently approximate. As to the 
proper decay model, the main problem is that it is essentially non-relativistic, which is 
of course in principle very far from the real situation. It is known from quite a long time 
that quite surprisingly, non-relativistic decay or emission models may work well, but their 
accuracy cannot be estimated a priori, it has always to be judged a posteriori. Decay 
models must be necessarily combined with potential models giving the wave functions 
that must be folded into their general structure. In view of the rather naive status of 
the ^Pq model, we do not find it appropriate to use a sophisticated set of wave functions, 
but rather a simple-minded one, as explained in section 3.2. We must underline however 
that the oscillator radii that are used are not at all free parameters: they have to be fixed 
on the actual spectrum. On the other hand, the model contains free phenomenological 
parameters, namely the mixing angle of Ki states, 6'^-^, and the quark-pair-creation con- 
stant 7. These parameters are to be adjusted on the strong decay experiments themselves 
(additional information on the mixing can be obtained from the mass spectrum, or other 
types of decays). 

4 The plan of the paper 

In Section 2, we present a brief summary of the present status of the experience concern- 
ing the i^i-mesons. In Section 3, after having discussed the basic question of the mixing, 
we introduce the formalism of the theoretical model, namely the ^Pq quark-pair-creation 
model, used to predict the partial wave amplitudes for the quasi-two-body decays of the 
f^i-meson. In Section 4, we establish the general relation between our model predictions 
and the most extensive experimental results obtained by ACCMOR collaboration, which 

^See [6] for another approach based on phenomenological Lagrangian. 
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use the K-maXiix formalism to analyse the partial waves. We describe some of the prob- 
lems we have observed, which include the definition of the total iiTi-width, the phase 
space and threshold effects, the strong phase between different intermediate resonance 
("isobar") states. In Section 5, combining the experimental results on the i^i-decays 
and the predictions of the quark-pair-creation model, we determine the phenomenological 
parameters of this decay model, the Ki mixing angle 6ki and the universal quark-pair- 
creation constant 7, and we present the resulting numerical predictions. We compare 
our model predictions and the measurements of the ACCMOR and Belle collaborations. 
We also discuss the issues of the relative strong "offset" phases and the controversial 
i^ri(1270) — )■ KIT channel. We give our conclusions and perspectives in Section 6. 

2 Overview of the previous experimental i^i-decay 
studies 

Here we summarise the experimental results of the axial vector i^'i-resonance study. 

1. Two close in mass axial- vector mesons, i^'i(1270) and fri(1400), were disentan- 
gled in the experiments on the diffractive production of the I^^Ktttt) system in 
the Kp — )■ KiTTTp reaction, first by the group at SLAC [7] and then by the AC- 
CMOR collaboration in WAS experiment at CERN [8]. They also observed sepa- 
rately: one i^i(1270) in the strangeness-exchange reaction n^p AKtctc [9] and 
the other Ki{14:00) in the charge-exchange reaction K~p K Tr+vr^n [10]. In 
our study we rely mainly on the diffractive reactions which allow a more detailed 
study. The relative ratios of two dominant channels, K*^ and Kp, indicate that 
7^1(1400) decouples from the Kp, while the Kp decay mode of i^'i(1270) is domi- 
nant (see Table 1). This decay pattern suggests that the observed mass eigenstates, 
ii'i(1270) and i^i(1400), are the mixtures of two strange axial-vector SU{3) octet 
states Kia(^Pi) and Kib(^Pi), as explained later. 

2. The i^'i-resonances were also observed and studied in r-decays, r — Kiv^, by 
TPC/Two-gamma [11], ALEPH [12], OPAL [13] and CLEO [14] collaborations. 

3. Radiative i?-decays involving the i^i-mesons were also observed by the Belle col- 
laboration [5]. The data indicate that B{B fsTi (1270)7) > ^{B ^ i^i (1400)7). 

4. Quite recently the Belle collaboration published a paper on B ^ J /il){il)')KTnT 
decays [15], which will be discussed in detail later. 

5. In addition, the BABAR collaboration reported the measurement of the branching 
ratios of neutral and charged i?-meson decays to final states containing a -ft'i(1270) 
and ii'i(1400)-meson and a charged pion: 5° iri(1270/1400)+7r" and B+ 



7 



i^i(1270/1400)°7r"'' [16]. In order to parametrize the signal component for the pro- 
duction of the i^i-resonances in 5-decays, the i^-matrix formahsm, used in the 
analysis by Daum et al. in [8], was applied for the model description. Since only 
some parameters, used in the analysis of the ACCMOR collaboration, have been 
reported, the BABAR collaboration refitted the ACCMOR data in order to deter- 
mine the parameters describing the diffractive production of the i^'i-mesons and 
their decays. One observes that some results are somewhat different. 





j^ACCMOR 

GeV/c2 


rACCMOR 

MeV/c2 


B{K*7r)s 




BiKp)s 




1.27±0.007 


90±8 


0.13±0.03 


0.07±0.006 


0.39±0.04 


Ki{iAm) 


1.41±0.025 


165±35 


0.87±0.05 


0.03±0.005 


0.05±0.04 



Table 1: Fitted masses, total widths and partial branching ratios of Ki{l^) decays into 
vector-pseudoscalar states, measured by the ACCMOR collaboration in the Kp Kmrp 
reaction for the low momentum transfer to the recoiling proton [8], and tabulated in PDG. 
The total widths seem to be misleading for the calculation of partial widths, as discussed 
later in the text. 



3 The theoretical model 

Before presenting the ^Pq, we begin by explaining the question of the mixing of Ki sates, 
which is a basic assumption of all the approaches, since there is no theoretical approach 
predicting quantitatively this mixing. 

3.1 The mixing of the Ki resonances 

In the quark model there are two possible states for the orbitally excited axial-vector 
mesons: J^*^ = 1^+ and J^'" = depending on different spin couplings of two 

constituent quarks. In the S'f/(3)-limit these states do not mix in general, but since 
the s-quark is actually heavier than the u- and rf-quarks, the observed i^i(1270)- and 
1400) -mesons are not pure l^Pi or l^Pi states. They are considered to be mixtures 
of non mass eigenstates Kia and Kib- Introducing a Kia — Kib mixing angle 6'xi, mass 
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eigenstates can be defined in the following way [17] 



17^1(1400)) = \Kia) cos9k, - \Kib) sinOK, 



(2) 



Since all of SU (3) operators can be expressed as combinations of isospin, U- and y-spin 
operators, if an operator describing the interaction is invariant under the S'?7(3)-group 
transformations, it is also invariant under the isospin, f/-spin and l^-spin transforma- 
tions [18]. However, it is sufficient to require the invariance only under the isospin and 
f/-spin (or \^-spin) transformations, since \^-spin is dependent on the isospin and f/-spin 
and the l^-spin operators can be obtained from the f/-spin operators by an isospin trans- 
formation ([/-spin can be turned into \^-spin via rotation by 120°). 



Analogously to G-parity, one can define U- and V^-parities: Gu = C{—1) and G 



V 



C(— 1)^ respectively, where C is the charge-conjugation parity of the neutral non-strange 
members of the multiplet. The neutral and charged kaons in the octets are the eigenstates 
oi U- and V^-parities and always have U oi V = 1 respectively. 

In the S'[/(3)-limit two kaons that belong to the octets of the same spin but opposite 
C-parity can not mix. To illustrate it, one can consider a matrix element of some arbitrary 
operator O between two neutral kaons from different octets [19,20]: 

{Ka\0\Kb) = (KaIG^j^GuOG^'GuIKb) = CACB{KA\GuOGij'\KB) (3) 

If the O operator is S'f/(3)-invariant, i.e. G^OG^^ = O, the matrix element of the 
transition {Ka\0\Kb) = unless Ca = Cb- 

Strong interactions can break the S'f/(3)-symmetry and produce the mass splittings. 
It is experimentally confirmed that isospin is conserved in strong interactions. Hence, 
if the strong interaction operator breaks the SU (3)-symmetry, U- and \^-parities are 
not conserved anymore, even if G-parity is conserved. In this case GuOG~[j^ ^ O and 
consequently {Ka\0\Kb) 7^ and the mixing takes place. 

That this mixing is indeed the effect of the symmetry breaking can be explicitly seen in 
quark models. At the level of bound states of a potential model. It is induced, for instance. 



•^To be able to compare with other mixing angle estimations, one has to be careful due to the 
different parametrizations that are used in the literature. For instance, in the analysis by Carnegie 
et al. [7] the parametrization is |ifi(1270)) = \Kia) cosO'^j^^^^^^ + \Kib) siii0'j^^^'^^\ |i^i(1400)) = 
— \Kia) sinO'"^^'^'^'' + \Kib) cos9^^^'^'^\ To compare with the results made by Damn et al. [8], 
parametrization is written as follows: |i^i(1270)) ~ — ji^i^i) sin 6'^^'^'^*^'^^'' + \Kib) cos9^j^^'^^^'~^^\ 
|i^i(1400)) = \Kia) cos f?^*"*"*^'^'^^ + \Kib) sin 6*^*" ^^'^■'^^ Comparing the fitted effective couplings one 
can see that the coupling to Kib has a different sign in these two definitions. Since one can measure only 
the absolute value of the amplitude, this sign changes nothing and hence it is possible to redefine the sign of 
this coupling in the paper by Daum et al.. After that one can easily establish the correspondence between 

(ACCMOR) (SLAC) 

these two forms of parametrization and the one we use in this paper: Bk^ — 9p^^ = 90° — 6]^^ 
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by spin-orbit forces with different s and u, d quark masses. A mixing is also generated 
by the two-meson loops, due to quark pair creation and annihilation in the bound states, 
as is explained below in the /T- matrix approach, subsection 4.1. In this approach, the 
real mixing must be understood as the one of the i^-matrix couplings corresponds to the 
effect of the real part of the loops, while additional, complex, mixing would be present in 
the physical couplings. 

One can see easily why the loops and the SU (3) breaking generate mixing. For in- 
stance, the K*^ and Kp loop contributions connect the Ka and the Kb, since both states 
are coupled to these channels. In this way it generates the mixing. The two contributions 
cancel each other if one sets Mk* = Mp and = m^, i.e. in the case of the exact 
S't/(3)-symmetry. It must be emphasized however that this mechanism of loops does not 
lead to the actual calculation of the mixing angle, because one would have to sum over a 
very large number of possible intermediate states. Therefore, in this approach, it remains 
an independent phenomenological parameter, which has to be fixed through confrontation 
with data. 

It is to be noted that, anyway, no fundamental calculation of the mixing has been 
produced. 

3.1.1 Previous phenomenological determinations of the mixing angle 

Here, we gather all the various estimations of Ok-i^. 

On the other hand, there have been, in the past, many attempts to determine the 
mixing angle, both from experimentalists and theoreticians, but in both cases only through 
phenomenological analyses. The phenomenological analyses have concerned the masses 
(with additional assumptions, since SU (3) alone does not enable to fix the mixing angle 
from the masses), the t ^ KiVj. decays, the B ^ Ki transitions, and, mainly, the strong 
decays Ki — )■ K'K'k through K*Ti and Kp channels. Indeed, the pattern of the latter is 
very sensitive to the mixing angle. 

The angles are given according to the definition above, Eq. (2). However, one must 
warn that it does not completely fix the definition, since there may be different choices of 
the phases of the states. In general, it is difficult to establish completely the connection 
with our own definition in the present paper, so we only state the absolute magnitude of 
the angle. 

• In the experiment, carried out at SLAG by Carnegie et al. [7], the mixing angle was 
determined from the SU{'i) couplings to the K*tt and Kp channels to be 9ki = 
(41 ± 4)°. On the other hand, the partial wave analysis of the WA3 experiment 
data, done by the ACCMOR collaboration (Daum et al. [8]), gives Ok^ = (64 ± 8)° 
and = (54 ±4)° for the low and high momentum transfer to the recoiling proton 
respectively. 
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In the reanalysis of ACCMOR data by BABAR [16], using the low t-data, the 
refitted value of the Ki mixing angle turns out to be 72° compared to 64° from the 
ACCMOR fit. 

In the work by Suzuki [17], the mixing angle is determined by three different ap- 
proaches. One is in order to explain the observed hierarchy in the Ki strong decays 
to K*^ and Kp, like has been done by SLAC and ACCMOR. Another is the S'?7(3) 
analysis of the masses of the two octets, but with additional assumptions . Finally, 
the suppression of r — )■ Ki{lAQQ)u with respect to r — )■ f('i(1270)z/ is considered. 
Two possible solutions for the Ki mixing angle were found: ~ 33° or 57°. 

In the analysis of r — )■ Kiv done by the CLEO collaboration [14], the Ki mix- 
ing angle was determined from the measured ratio B{t — Ki{12'1{))Vt-) / B{t — )■ 
Ki{lAm)Ur)- 0Kr = (69 ± 16 ± 19)° for 5 = 0.18 and Ok, = (49 ± 16 ± 19)° for 
6 = —0.18 where \6\ = "j-"-^^ ~ o.l8 is a phenomenological SU{3) breaking 
parameter. This result is consistent with the calculation by Suzuki [17]. 

In the work of Blundell, Godfrey and Phelps [21] 

— The mixing is discussed using the results of the TPC/Two-gamma collabo- 
ration: B{t- fs:i(1270)-z/,) = (0.41^°:^^)% and B{t- Ki{U00yUr) = 
(0.76]l^g33)%. This would seem to mean that the rate into ifi(1400) is larger 
than into A'i(1270), although their errors are too large to make a strong state- 
ment. Anyway these numbers have been superseded by the CLEO data, which 
show the contrary. 

— The strong decays of the /^i-mesons to the final states K*n and Kp were 
studied as well in order to determine the mixing angle. A fit of the ex- 
perimental data on the partial decay widths F (7^^1(1270/1400) — )■ K*tt) and 
F(i^'i(1270/1400) — )■ Kp) was used for the 6'7^^-determination. 

* Performing a x^"fit with the predicted decay widths, calculated within 
the pseudo-scalar-meson-emission model, using simple harmonic oscillator 
wave functions with a single parameter /3 = 0.40 GeV, the fitted value of 
the mixing angle was obtained to be O^i = (48 ±5)°. 

* The strong i^i-decays were also calculated using both the fiux-tube-breaking 
model and the ^Pq model for several sets of meson wave functions. In all 
cases a second fit was performed by allowing both 9k-^ and the quark- 
pair-creation constant 7 to vary, what reduces the significantly. Using 
simple harmonic oscillator wave functions with (3 = 0.40 GeV, comparison 
of the predicted decay widths by the ^Pq model to experimental results 
gives 6ki = (45 ± 4)°, while the fiux-tube-breaking model's prediction 
gives = (44 ±4)°, both appreciably different from our central value 
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— 60° with the same set of wave functions. Their last resuh for 
is shghtly changed for the case of use of different set of the meson wave 
functions from Ref. [22]: 9k, = (51 ± 3)°. 

• A detailed study of the B A"i (1270)7 and B -> 7^1(1400)7 decays in the light- 
cone QCD sum rules approach was presented by Hatanaka and Yang in [23]. The 
sign ambiguity of the mixing angle is resolved by defining the signs of the decay 
constants fx,^ and f^^^. 

— From the comparison of the theoretical calculation and the data for decays 
B — 7- Ki'-f and r — > Kiu^, it was found that 6ki = —(34 ±13)° is favoured 
within the conventions of Hatanaka and Yang. It is difficult to establish the 
relation with our own convention as regards sign. 

- The predicted branching ratios, B{B JsTi (1400)7) and B{B ^ i^i (1270)7), 
are then in agreement with the Belle collaboration measurement within the 
errors. 



In summary, the cleanest way to extract the mixing angle is certainly, in principle, the 
determination from the ratio of B{t — > i^i(1400/1270)z/^), if the data were sufficiently 
accurate. At present, we believe that the best way remains the study of strong decays, as 
we do in this paper. 

3.2 The Quark-Pair-Creation Model 

There are several additive quark models of strong vertices. All these models relate to 
the recoupling coefficients of unitary spin, quark spin and the quark orbital angular mo- 
menta, but differ in the dynamical description. One of the simplest additive quark model 
describing three- meson vertices is the naive quark-pair-creation model (QPCM), with 
a ^Pq structure for the pair, formulated by Le Yaouanc, Oliver, Pene and Raynal [24] 
starting from ideas of Micu and of Carlitz and Kislinger [25,26]. The model has then 
been extensively applied and discussed by many authors, including the same authors 
(see Ref. [27] and some references therein) and the group of N. Isgur in Canada (for in- 
stance Refs. [28,29]). As in the usual additive quark models with spectator quarks, the 
quark-antiquark pair is "naively" created not from the ingoing quark lines but within 
the hadronic vacuum. The strong interactions vertices in the QPCM are expressed in 
terms of the explicit harmonic oscillator spacial SU{Q) wave functions (compared to the 
work by Micu [25], who just fitted the various spacial integrals using the measured decay 
widths, what does not allow to study the polarization effects) and a nonlocal vacuum 
quark-aniquark pair production matrix element, depending on the internal quark mo- 
menta (while Carlitz and Kislinger [26] neglected the internal momentum distributions). 
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Contrary to the QPCM by Colglazier and Rosner [30], the ^Pq structure of the created 
pair describes any decay process of any hadron, using one universal parameter. The other 
model parameters are those of the hadrons themselves (potential model), and not relative 
to the decay process as in [30], where the various extra couplings between the pair and 
the incoming meson depend on the nature of the hadron states and may be weighted by 
different arbitrary coefficients for different hadrons. 

The naive QPCM has the advantage of making definite predictions for all hadronic 
vertices and moreover, contrary to the other works, it predicts the relative signs of the 
couplings. Another appealing feature of the model is that it consists only one phenomeno- 
logical parameter (the quark-pair-creation constant), what allows a much more general 
description and relates the amplitudes of different processes. The main weakness of the 
QPCM is that the emitted hadrons are considered to be non-relativistic. Thus one has to 
look for the decays that are not significantly sensitive to these effects. 

A specific study of the strange axial- vector mesons was first done by Blundell, Godfrey 
and Phelps [21], who studied the properties of Ki by combining wave functions inspired by 
the Godfrey-Isgur quark model to describe the bound states, and the fiux-tube-breaking 
or ^Pq models to describe the decays. Although we start from the same basic ^Pq model, 
we give a much more extended study, which is, in particular, required for the purpose of 
the A-^ determination. We make a rather different discussion, especially, for the relation 
between theory and experiment. We clarify the relation with the ii'-matrix analysis, which 
is the tool used by the main experiment, that is the ACCMOR experiment. We discuss 
the definition of widths, which appears very ambiguous due to threshold effects. We also 
include a treatment of the off-shellness (i.e. damping factor). In addition, we explore the 
system of phases, which is one main achievement of the model (as well as it has been in 
the baryon decays). Finally, we discuss in detail the most problematic ktt channel. These 
differences will become apparent from the rest of the paper. 

3.2.1 Formalism 

In the QPCM, instead of being produced from the gluon emission, the quark-antiquark 
pair qq (see Fig. 1) is created anywhere within the hadronic vacuum by an operator 
proportional to (uu + dd+ss) S-p where S refers to spin 1 and p is the relative momentum 
of the pair. It is combined with the initial quark-antiquark system g2^i and produces 
the final state B{qiq)C{qq2)- The initial spectator quarks are supposed not to change 
their SU{3) quantum numbers, nor their momentum and spin. In order to conserve the 
vacuum quantum numbers the pair must be created in the ^Pq state due to P = — (— 1)^ 
and C = (—1)^"'''^ parity conservation with 0-total momentum (^3 -|- ^4 = 0) and to be 
a S'f/(3)-singlet. Thus the matrix element of the quark-antiquark pair production from 
the vacuum is unambiguously constructed with the help of the spins and momenta of the 
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quark and antiquark only [24]: 



(gg|f,„,|0) = 5ih + hh 5^(1, rn; 1, -m|0, 0)yrih - h)xI"'<Po (4) 

m 

where 7 is a phenomenological dimensionless pair-creation constant (which is determined 
from the measured partial decay widths and taken to be of the order of 3-5), Xi"^ ^ire the 
spin-triplet wave functions, 0o = -^{uu + dd + ss) is the S'f/(3)-singlet and 3^™ represents 
the L = 1 angular momentum of the pair. 

Taking the matrix element of the pair-creation operator between the SU (6) harmonic- 
oscillator wave functions of hadrons, the matrix element for the decay A ^ B + C can 
be written as: 

{BC\f\A) = 7 ^(1, m; 1, -m|0, O)*^*^^'^';^^^''''^ (5) 

where $ = xT^ SU{Q) spin-flavour wave functions and Im^'"^ are the spacial in- 

tegrals dependent on the momentum of the final states, which are computed in Appendix. 




Figure 1: Three-meson vertex in the quark-pair-creation model. 

Assuming A, B and C to be an axial vector, pseudoscalar and vector mesons respec- 
tively, the spin part of the matrix element can be written as 



XcXBXAXpair = ^(^,mi; ^,m3|0,0)(^,m4; ^,m2|l, Ac) 

(6) 

X (2'"^i5 ^,m2\SA,ms^)il,TnL^; SA,ms^\l, XA)i^,rn4; -,m3|l,-m) 

Consider for instance ii"*°7r"'" decay mode of i^i-meson. After the summation over the 
spin projections the calculated helicity amphtudes for the Kia (1^-Pi) and Kib (1^-Pi) 
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will be (the definition of the helicity amplitudes and their relation with the partial wave 
amplitudes can be found in Appendix): 



•'^loo ~ 6 ' -^10 ~ I Q 

The corresponding amplitudes for the K~^p^ mode are obtained by multiplying the K*^7r^ 
amplitudes by l/\/2 and changing the sign of i^'iA-part. 

Taking into account the isospin factors for different charge states \ the generalized 
amplitudes are summarized in Table 2. The functions S and D are defined as 

PTn AABC) AABC) AABC) AABC) 

W 2 18 ' ' V 2 18 ^ ' 



Decay mode 


As 


Ad 






-V2D(-^i-^*-) 


KiA ^ K*n 




_JJ{KiK''it) 


K,B ^ Kp 






KiA ^ Kp 


v/2^(^i^^) 





Table 2: Partial wave amplitudes of Kia{1^Pi) and Kib{1^Pi) decays into vector- 
pseudoscalar states, calculated within QPCM. 

One has to point out that our treatment obeys the SU (3)-symmetry. SU{'i) breaking 
effects are present only in two places: 1) we use the physical observed masses of hadrons 
to calculate the momentum transfer of the decay and the phase space; 2) we introduce 
mixing between the Kia and Kib states. 

Then the decay amplitudes of the physical Ki states into K*tt or Kp final states can 
be expressed as functions of the pseudoscalar meson momentum in the Ki reference frame 

''The amplitudes were calculated for K*^Tr~^ and — >■ K~^p^. The amplitude of Kp must 

be divided over \/2 due to isospin wave function of To obtain the general amplitude which doesn't 
depend on the charge combination one has to divide over the isopin factor: —^/2/3 for K* and -^/l/S for 
p since for the matching with the relativistic form factors the charge combination is not relevant. Finally 
one obtains the factor \/3/2 in Eq. (8). 



15 



and the mixing angle 6ki- 



As{Ki{1270) K*7T/Kp) = S^^'^*''^^'^p\V2smeK, T cosOk,) 
Ad{Ki{1270) K*Tx/Kp) = D'-^^^*''I^^^p\- s^Ok, T v^cos^/^, 



(9) 

As{Ki{lAm) ^ K*7T/Kp) = S^^^^'^'^^^f\V2coseK^ ± sin^^^J 



AniKiilAm) ^ K*Tx/Kp) = D''^^^*''I^'^p\- cosOk, ± v^sin^^J 

Correspondingly, the partial decay widths can be determined by using amplitudes squared 
from the Eqs. (9) multiplied by the phase space factors: 



r?;,^-(ir, VP) = 8-'^^\AsMK. -> vp)\\ (10) 



Note that all the signs in the expressions for amplitudes have sense only within defi- 
nite specific conventions. The ones in our work are defined in Appendix B. On the other 
hand, the signs of the products of the couplings of the two successive decay processes 
from the same Ki state, that is, when we multiply by the decay amplitude of the iso- 
bar, make sense and the relative signs are observable because the final state Kinr is the 
same, and all phase arbitrariness cancels. It is an important feature of the model that 
it can predict all these observable signs. As will be seen in subsection 5.3, these predic- 
tions are remarkably verified by experimental data. 



3.2.2 The choice of the spatial wave functions 

The unknown parameters of the model are the quark-pair-creation constant 7 and the 
Ki mixing angle, which we determine by fitting the experimental data on the i^'i-decays 
(see the next section). However, before proceeding to this determination, the model must 
be specified by the choice of the set of meson wave functions. In accordance with a 
fact that the ^Pq model is a simple model, we will remain within the traditional SU{6) 
approximation which describes rather well ordinary radiative decays (e.g. u — )■ 717). This 
includes the 5't/(3)-symmetry approximation which anyway is also present in the ^Pq 
model through the fact that the quark-pair-creation constant is the same for all reactions. 
In this approach the effect of the SU{3) breaking is taken into account only through 
the dependence of the decay momentum of the physical hadronic masses. For practical 
reasons, we choose a set of harmonic oscillator wave functions, which are known to give a 
reasonable approximation. 

Here one has to stress that the harmonic oscillator radius of the meson wave func- 
tion {ip{r) oc exp(— r^/2i?^), for details see Appendix B) is not a free phenomenological 
parameter. In principle, it can be predicted by the quark-potential model describing the 
bound states of two quarks. To get a first and rough estimate we can use the following 
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relation, obtained in the non-relativistic harmonic oscillator model for the energy shift 
between the ground state and the first radial excitation: 



with nig being the quark mass, which can be standardly estimated from the magnetic 
moment of the proton: fip = = Whence ~ 0.34 GeV. AEi can be estimated 
from the energy of the L = 1 state of the order of (1.2-1.3) GeV and the weighted average 
energy of the ground state (3mp + m7r)/4 ~ 0.6 GeV. Then the estimated radius is given 
by 



R=J- Wt ^ ~3GeV-^ (12) 

y A^im, y (1.25 - 0.6) 0.34 ^ ' 

On the other hand, it is obvious that this approximation of the Schroedinger equation 
with the harmonic oscillator potential is rather naive: the realistic potential is known to be 
of the form of linear (that describes confinement) plus Coulomb potential. One has also to 
notice that the application of the use of the non-relativistic character of the Schroedinger 
equation to the heavy-light systems is dubious. Therefore, one could take a value inspired 
by the well known model of Godfrey and Isgur. Of course, in the latter model the solutions 
are no longer the harmonic oscillator wave functions. However, such harmonic oscillator 
wave functions can represent a good approximation if the radius R is adjusted. For most 
L = 0, 1 states one finds in this model the typical value R ~ 2.5 GeV^^ [28]. For our 
predictions we therefore adopt a set of wave functions with a common harmonic oscillator 
radius having precisely this value, 

R = 2.5 GeV"^ (13) 

This is one of the choices made by Blundell et al. [21]. We must warn that in the model 
of Godfrey and Isgur, pion and kaon have actually quite smaller radius (~ 1.4 GeV~^ 
[28]) due to the strong spin-spin interaction force. If we were adopting the low values 
for the Goldstone boson we would obtain unsatisfactory results. For example, using 
Rt, ~ 1.4 GeV~^, we can not reproduce correctly the D/S ratio in the hi — )• uj-k decay 
which is precisely measured {D/S = 0.28). The use of the exact wave functions of 
the model of Godfrey and Isgur [22] does not seem to improve the situation; one finds 
D/S = 2.5/14 ~ 0.18 from the tables of Kokoski and Isgur [28]. 

Of course, although it has not been commented in previous works, this fact is disturb- 
ing, since the spin-spin force is present indeed in spectrocopy, and therefore, it should be 
more realistic to include its effect. In addition to empirical success, the choice of equal 
radii can be motivated in the spirit of the SU{6) approach. It must be remembered in- 
deed that old quark model very naive calculations have succeeded well with this SU{6) 
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symmetry, for instance to relate a; — )■ 717 to magnetic moments. Now, the ^Pq model is 
also in this very naive spirit: it is non-relativistic in essence. 

3.3 The issue of the damping factor 

In the end of this introduction of the theoretical model, we discuss the necessity of in- 
troducing an additional cut-off on momenta (or damping factor) in the coupling vertices. 
Generally speaking, there is need of a strong cut-off for calculations involving far off-shell 
particles, once the model has been adjusted on real decays. Indeed, the natural fall-off 
provided by simple continuation of the ^Pq model, due to the wave functions, is seen to 
be much too weak. The need for this cutoff appears in various circumstances: 

• In the branching ratios, obtained by the integration over a large phase space, like 
for the production of Kt^tt (e.g. B — )■ Rtttt'^, ip) or similar. For instance. Belle [15] 
defines branching ratios by the ratios of integrals over the whole phase space. If there 
were not such a cutoff, the higher partial wave contribution like D-waves would be 
found much too large with respect to S waves, due to the centrifugal barrier factors 
fc^^, which increase too much at large mass of the Knn system ^. 

• Departure of the resonance line shape from the Breit-Wigner formula. Resonances 
are usually described by multiplying the standard Breit-Wigner (and the width) 
by the so-called "centrifugal barrier" factors. The term is ambiguous, since these 
factors includes both the proper centrifugal barrier effect, which is the universal k"^^ 
automatically present in partial waves (increasing with the momentum), and the 
damping factor, which is highly model dependent, and decreases with the momen- 
tum. In fact the prototype of such factors are the Blatt-Weisskopf factors of nuclear 
physics, also commonly used by experimentalists in particle physics. They are de- 
duced for a spherical well potential, which is obviously very naive. One consequence 
of this particular set of factors is that there would be no damping for S waves, which 
is not true in more realistic models (harmonic oscillator wave functions in the ^Pq 
model give a Gaussian damping in all waves). 

The accurate studies of the resonance shapes show directly a departure from the 
standard Breit-Wigner shape, e.g. for the A(1236) [31] or the K*(890), see Ref. [10]). 

• Contribution of loops to the self-energy 

The need for the cutoff is also shown by calculations of the hadronic loop contribu- 
tion to the self-energy of mesons (see subsection 4.1.2), which involves integration 

^Experimentally, the problem does not appear in the work of Belle, because they do not introduce 
D-waves for the iiTi-decays. 
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over the possible momentum up to infinity. In the ^Pq niodel [32], the contribu- 
tion to the self-energy would be much too large for D-waves, in spite of the cutoff 
naturally provided by the gaussian wave functions, yielding finally a bad spectrum. 

One obtains a natural damping factor through the Gaussian factors e~^^'^: 

oc (3 - ak'')e~^''\ Ad oc ak'^e-^''' (14) 

but one finds /3 ~ 0.3 GeV^^ which is much too small ; it does not reduce effi- 
ciently the D waves contributions for the loops and neither for the off-shell situa- 
tions we consider. Following Ref. [32], we introduce the empirical Gaussian cutoff 
exp[— /3'(A;^ — k^)] with f3' ^ 3 GeV^^, where ko is the decay momentum when all 
the particles are put on-shell: 

As oc {3- ae)e-^'' x e-P'^''~'^\ An oc ak\-^^' x e~^'^^'-^'-^ (15) 

With this additional damping factor one finds that the integrated D / S-ratio becomes 
stable. The isobar {K* / p) decay does not depend much on the damping factor. However, 
another effect then appears in the decay rate from the parent Ki to one isobar and 
one stable particle: integrating over the mass of the isobar, the calculated partial width 
depends on the presence of the damping factor for the decay of the parent Ki resonance 
to an off-shell isobar. The low end of the isobar mass spectrum corresponds indeed to 
large off-shell momenta. This effect has been duely taken into account in our calculations. 

In the calculation of presented in our previous paper [3] . the effect of the introduc- 
tion of this damping factor in the decay amplitude of the Ki is important. Indeed, the 
interference of several channels needed to obtain a non-zero imaginary part of n- {J x J*) 
requires a large off-shellness of the intermediate isobars. We find that this quantity is sen- 
sitive to the presence of the D waves, and then to the introduction of the damping factor. 

4 How to compare the theoretical model computa- 
tion with the experimental data? 

Let us stress that the use of experimental data in our work is twofold: first determine the 
model parameters 7 and ^i^^, and then check the validity of our model predictions. 

In this section we will explain how one can relate the quark model predictions for the 
decay partial widths of Ki to the /^-matrix analysis of Daum et ai, which is the main 
source of experimental information. 

Indeed, the main experiments on the i^i-decays [7, 8] were analysed with the same 
i^T-matrix formalism developed by Bowler et al. [33] and obtained very similar results. 
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We use in our analysis the parameters of the analysis done by Daum et a/.(ACCMOR 
experiment) which seems to be the most detailed. On the other hand, there are certain 
physical parameters of the fit which are not tabulated in the this paper. Then we also 
use, where necessary, the results of the iiT-matrix re-analysis of the ACCMOR data by 
the BABAR collaboration [16]. 

Let us now emphasize that the very extensive work of Daum et al. consists of two 
distinct steps: 

• The first one is the partial wave analysis (PWA) where the Ktttt three-body final 
state is decomposed into a sum of quasi-two-body "partial waves" {K*7i, Kp, etc.) 
with various quantum numbers of the total spin and orbital momentum. In this first 
step there is no reference to any parent resonance like Ki. This step corresponds 
to the fitted values of the quasi-two-body partial wave amplitudes plotted with the 
corresponding error bars in [8] . 

• The second step is the fit of the partial wave amplitudes, extracted on the previous 
step, within the i^-matrix formalism in order to study the structure of the initial 
parent Ki resonance and its properties (pole masses, couplings to various decay 
channels, etc.). 

Let us stress that this two-step procedure is different from the modern Dalitz plot 
analyses where the isobar and parent resonances are included together in one unique 
formula of the total amplitude. In that case the total amplitude is written directly as a 
product of the parent resonance decay amplitude and the amplitude of the subsequent 
decay of the isobar taking into account the width effects of the unstable resonances by 
the Breit-Wigner forms. 

We do not question the first step; we rather indicate various difficulties which we have 
encountered in trying to use the /^-matrix parameters from the analysis of Daum et al.. 
In the following subsection, we first recall the general i^T-matrix formalism and then its 
relation to the quark model. 



4.1 The ii'-matrix formalism and the quark model 

In order to extract our theoretical parameters, 7 and 6ki, we need the experimental partial 
widths. We also need them to verify our prediction of the model. And the question is: how 
to define a partial width? Resonances are often parametrized in terms of the Breit-Wigner 
form 

BWl:^^\m) oc ^, or BWj:^\m) oc 5^- — (16) 

in the non-relativistic and relativistic cases respectively. Resonance width, in principle, 
depends on energy, Tr{m). This approximation assumes an isolated resonance with a 
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single measured decay. If there is more than one resonance in the same partial wave 
which strongly overlap, an elegant way that provides the unitarity of the S'-matrix is to 
use the i^-matrix formalism for the two-body decays of the resonance states (for more 
details see Appendix) ^. 

4.1.1 General definitions in the _ft'-matrix formalism 

From the unitarity of the S'-matrix 

S =l + 2ip^2Tp'2 (17) 

one gets 

T -T^ = 2iT^pT = 2iTpT^ (18) 

where the diagonal matrix Pij{m) is the phase space factor which is discussed in detail 
later in this section. In terms of the inverse operators Eq. (18) can be rewritten as 

(T^)-^ -T-^ = 2ip (19) 
One can further transform this expression into 

(T-i + ipy = T-^ + ip (20) 
Using the definition of the i^-matrix 

K-^^T-^^ip (21) 
one can easily find from Eq. (20), (21) that the ii"-operator is Hermitian, i.e. 

K = K^ (22) 

From the time reversal invariance of 5* and T it follows that K must be symmetric, i.e. 
the K-matrix can he chosen to he real and symmetric. Resonances should appear as a 
sum of poles in the i^-matrix. In the approximation of resonance dominance one gets 
therefore 

^ y faHfa'j ^23) 

nia' —m 

a' 

where the sum on a' goes over the number of poles with masses m^/. In the common 
approximation in the resonance theory, the couplings are taken to be real. 



^Note that in the case of two overlapping resonances the Breit-Wigner parametrization of the amph- 
tude satisfies the unitarity condition of the S'-matrix only with the complex couplings satisfying certain 
condition. As we demonstrate later, these complex couplings can be obtained from the real if-matrix 
couplings by a complex rotation. 
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The partial and total i^-matrix widths can be defined as 



(24a) 
(24b) 



Note that the i^'-matrix width does not need to be identical with the width which is 
observed in experiment nor with the width of the T-matrix pole in the complex energy 
plane. 

4.1.2 Relation between the couplings in the ii'-matrix formalism and the 
quark model 

In this section we identify in a systematic approach the couplings deduced from the ^Pq 
quark model, including the mixing of Ki resonances, with the couplings introduced in 
the i^-matrix formalism by Bowler et al. [33]. To justify this identification, we establish 
the connection between the formalism, introduced in the previous section, and the quark 
model. 

1. To make explicit the discussion in Ref. [34], we distinguish two types of interactions: 



• The first type of interactions is described by Hamiltonian Hq, which describes 
the qq potential of the bound states of mesons, {a°, 6°, . . . }. It generates the 
initial meson masses and wave functions which are used to calculate the matrix 
elements of meson decays in the quark model (see next item). 

• The second type of interactions, described by Hamiltonian H', represents the 
interaction vertices connecting these bound states to the continuum of all pos- 
sible states of two interacting mesons, {i,j, . . . }: 



We commonly call these vertex interactions "couplings" . These couplings can 
be precisely calculated within the ^Pq quark-pair-creation model. With ade- 
quate choice of phases of the wave functions of the bound states the couplings 
can be set to be real. 

2. No direct interaction is assumed between two mesons. Nevertheless, there is rescat- 
tering since a meson pair can annihilate into one bound state and then be created 
again from the decay of this bound state. This rescattering process can be iterated 
arbitrary number of times, what is equivalent to a resummation of meson loops 
between the initial and final vertices (see Fig. 2). 



(25) 
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Figure 2: Rescattering process. 

All these possible processes can be resummed into a matrix propagator HaOto con- 
necting two vertices. Let us call the scattering energy m. The couplings can be in 
principle energy dependent, i.e. be a function of m. This is the case with centrifugal 
barrier or damping factors, which are indeed present in our model. However, for 
simplicity of the presentation we assume them to be constant. Then, defining the 
"bare" scattering amplitude for the first diagram in Fig. 2 

m^o — m 

and resumming all possible digrams with leads to the scattering amplitude 



aO aO,k,W aO,k,bO,l,cO 

= ^ fiaO^aObofbOj 
aO,bO 

where the propagator defined as 



(27) 



(28) 



n<jOfeo(m) = h 7 TT r/fc + . . . 

TTLaO — m [mao — 'm)[mi,o — m) 
T(m) 

M{m) - - m 

with Ik being the loop integral for the rescattering loop of the k^^ channel. In the 
case where the couplings / depend actually on m, one should include the coupling 
factors relative to the internal lines of Fig. 2 in the loop integral. However, since 
we do not attempt to calculate actually the loops, we see that there is no need to 
introduce this complication. 

The mass matrix in (28) is: 

( M - ^ J = m^O JaOfeO - ^ fiaO fboJi (29) 
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It is in general non- diagonal. It contains 

• the initial diagonal mass matrix diag(maO, m^o, . . . ) of the bound states; 

• the contribution of the loops for each possible channel, which can be non- 
diagonal since common two-body channels can couple to two different bound 
states. The loop integrals contain real and imaginary parts, which appear only 
when a two-body channel is open at the energy m. 

3. Now the mass matrix must be diagonalized in two steps as explained in [34]. One 
first diagonalises the real part, M, then one passes to a diagonalization of the full 
new matrix, M' — iV /2. 

(a) Diagonalization of the real part of the denominator of n^OfjO matrix, M, leads 
to the introduction of the new diagonal mass matrix 



This mass diagonalization implies a simultaneous rotation of the couplings 
{fa°i}) leading to the new couplings {fa'i}- One thus passes to the masses 
and couplings of the ii'-matrix, Eq. (23). Of course, if there exists only one 
resonance which couples to the initial and final states, no rotation is needed. 
In this case all bare couplings {/aOj} coincide with the ones of the i^-matrix, 
{fa'i}- Thus, one can relate them with couplings calculated in the quark model. 
Otherwise, when there are two possible overlapping resonances, namely the 
two -R'l's, we have to make a rotation and introduce a mixing angle. We 
notice then that we have introduced an arbitrary rotation angle 9ki in our 
model computations which allows us to identify the set of the observed K- 
matrix couplings with the theoretical ones by the fit of data with our model 
predictions. This identification means that: 

• the effect of the real part of the loops, i.e. Re{Ik) in Eq. (27), are taken 
into account in our model; 

• mixing angle 9ki is not predicted by the model but is simply adjusted to 
data; 

• introduction of the mixing angle 6ki can also take into account the uncal- 
culated rotation of the pure spin states Kia and Kib into the eigenstates 
of Hamiltonian Hq due to the spin-orbit forces [35] . 

(b) The second step consists the diagonalization of the new mass matrix 




(30) 




(31) 
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This leads to the physical mass eigenstates and to the Breit-Wigner parametriza- 
tion with energy-dependent width. This new rotation that accomplishes the last 
transformation into the physical states must have a complex and the angle of this 
rotation must have a complex phase. This would lead to the complex couplings of 
the mass eigenstates to set of continuum states. As we have already mentioned in 
the text, this rotation seems to be rather small. 

4. Let us now discuss the dependence of various variables on the energy m. In principle, 
all the masses and couplings, produced by two previous steps are dependent on m 
because of the loop effects (the bare couplings themselves may depend on m, as 
is the case of our quark model, when one calculates the decay momenta, then the 
widths at the mass m. This then modifies the expression of the loop integral). This 
also implies that the real mixing angle 6ki is also energy dependent in principle. 
However, as regards the mass matrix, its real and imaginary parts have rather 
different behaviour depending on m. In first approximation, the real part of the 
mass matrix, which includes the sum of the large number of loops, varies slowly 
with m and can be considered as constants on a limited range of energy. This is 
what was done in the analysis of Daum et ai. On the contrary, the imaginary part, 
which corresponds to the partial widths of the opened channels, is a rapidly varying 
function near the threshold. 

One can go beyond the approximation of the real part of the mass matrix by taking 
into account that there is some variation near the threshold. This is obtained 
by analytic continuation of the phase space through the threshold. This effect is 
consistently included in the prescription of Nauenberg and Pais of the complex 
phase space, although we differ on other assumptions they made. This corresponds 
to having imaginary part of the widths generating a m-dependent mass shift. For 
instance, for the i^'-matrix width one have 



where the phase-space factor Pij{m) can be complex in general. 

Finally, one obtains for the physical states that the physical masses of i^'i (1270/1400) 
are varying slowly as functions of m while the physical widths are rapidly changing 
functions; moreover the mass of /ri(1270) has a more rapid variation around the 
peak due to the closeness of the i^'i-mass to the Kp threshold (see Fig. 10). 

In summary, one should identify the i^'-matrix couplings with the ones predicted in the 
^Pq model, with the real mixing effect included to define the initial states in this model. 
To establish the quantitative relation between the definitions in these two formalisms, we 




(32) 



25 



identify more exactly the partial widths, Fj^^- (Mpeak)'- 

r^ff ^'(Mpeafc) = r^,, = 2/i^,i?e[p,,(Mpeafc)] (33) 

where r^^j'"*^(Mpeafc) is the ^Pq model partial width, Eq. (10) ; F^^j, pij and fxii are 
respectively the partial width, the phase space (we use the real part of the phase space 
since pij is defined as a complex quantity as will be explained later) and the i^-matrix 
couplings in the formalism of Daum et al. Note that these are not exactly the common 
partial widths related to the Breit-Wigner analysis: the latter would be obtained by 
applying the complex rotation (see Eq. (31)). 

Now, Eq. (10) is valid only for the narrow isobar. If we have to take into account 
the effect of the finite width of the isobar, we have to integrate the quasi-two-body phase 
space over the Breit-Wigner of the isobar. One has to underline, that in this approach we 
do not have to integrate over the Breit-Wigner of the i^Ti-resonance unlike what is done, 
for instance, in Ref. [28]. We indeed calculate the width at the peak. On the contrary, 
if we would like to compare with the results of the Belle collaboration analysis [15], this 
approach must be changed and we would have to integrate over the whole three-body 
phase space of 5 — )■ KttttiIj to obtain the branching ratios. But, even in that case, it 
does not have sense, in our opinion, to integrate the decay widths themselves over the 
Breit-Wigner of the Ki. 



4.2 Observed problems in the experimental iC-matrix analysis 

As announced, we found several problems in using the experimental analysis: 

• Absence of the Kmr non-resonant contribution in the fC-matrix. 

We note that the A'-matrix of Daum et al. is composed only of two resonance poles. 
There is no non-resonant contribution which is usually parametrized as polynomial 
in terms of m in the AT-matrix parametrization. This implies the strong assumption 
that the quasi-two-body scattering of vector-scalar mesons {K*'7i and Kp) passes 
only through the Ki resonant intermediate states. 

• D-wave amplitudes issue. 

The results of the ACCMOR analysis show that the D-wave in A'i(1270) -> K*7r 
depends strongly on the production transfer t in the Kp — )■ Kirirp reaction. This 
fact may escape the attention of PDG reader, because it averages between two sets 
of data (low t, high t). As for the D-wave amplitude in the Kp channel, there is no 
information; only branching ratios are quoted in the paper but not the A'-matrix 
couplings and their phases which are crucial for our study. 
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• The problem of definition of the total width with threshold eff'ect 

When the mass of the resonance at the peak is close to a decay threshold, different 
definitions of the resonance width are no longer equivalent. Such possible definitions 
are the width at the peak r[Mpeak), the width at the S'-matrix pole, and finally the 
full width if measured at one-half the maximum height (FWHH) of the Breit-Wigner 
distribution defined as 

r™^m2-mi, (34) 
where mi and m2 are defined as two solutions in m of the equation 

fa'ib')lPlli^) _ 1 fa'(b')lPn{Mpeak) 

rria'ib') - m- iVa'(h'){rn) 2 ma>(h') - M^^ak - i'^a'(b')iMpeak) 

using the i^*7r channel (labelled as channel 1). 

The last two widths are found to be smaller than the first one. That is why the 
iri(1270) width, ri^,(i27o) = (90 ± 8) MeV/c^ [8], which is assumed to be defined 
as the full width if measured at one-half the maximum height of the Breit-Wigner 
distribution of Ki, is less by a factor 1.5-2 than the total width at the peak (see 
Table 3) which is computed using the i^-matrix couplings and summing over all 
possible intermediate channels, i.e. 



YP^f ^ 2J2 fhRe[p.^iMp^ak)] (36) 

i 

We find, indeed, for the later to be of the order of 200 MeV/c^ with the inclusion of 
the KTT channel (see Table 3). As a consequence, one observes a large discrepancy 
between the two possible definitions of the partial width that can be extracted from 
data of the ACCMOR collaboration: the partial width, defined in a "standard" 
way as r(/s:i(1270) ^ Kp) = Tk, x B{Ki{1270) ^ Kp), is less by a factor 2-3 
compared to the partial width at the peak, defined from the i^'-matrix couplings 
(see Table 4). The total width, defined by ACCMOR collaboration and tabulated 
in PDG, seems therefore to be misleading. It should not be used to compare with 
the quark model predictions. According to us, previous theoretical analyses (for 
instance, in Ref. [21]) unduely used for experimental partial widths the product of 
branching ratios with this total width of -ft'i(1270) quoted by PDG. 

• The problem of the phase space, pij 

In the expression of the T-matrix in the i^-matrix formalism the phase space factor 



27 





YACCmor^ MeV/c2 


MeV/c2 


YFWHh^ MeV/c2 


Ki(1270) 


90±8 


~190 


~80 


iri(1400) 


165±35 


-230 


-230 



Table 3: Experimental total decay widths, calculated using the fitted parameters from 
Ref. [8]. In our opinion, only the widths calculated at the peak must be used to compute 
partial widths from the branching ratios. Note that the D- waves are not included in the 
1 estimation. 



Decay channel i 


J- Kii — iJKii X J- i^Ti 5 

MeV/c2 


MeV/c2 


fs:i(1270) ^ {K*7r)s 


12±3 


28±26 


7^1(1270) ^ {Kp)s 


41±10 


122±28 


fs:i(1400) ^ {K*7r)s 


162±13 


211±59 


K,{UOO) ^ {Kp)s 


2±2 


20±25 



Table 4: Experimental partial decay widths, calculated using the fitted parameters from 
Ref. [8]. As is is underlined before, only the values from the last column must be used. 



Pij is defined as 

2ki{m) 

Pij{m) = dij (37) 

m 

Naively, ki, is the break-up momentum for the two-body decay channel i. But, in 
fact. Bowler et al. used for ki a particular formulation, proposed by Nauenberg and 
Pais [36], which tries to take into account two important effects: 

— The requirement of the analiticity of the amplitude. The simplest way to satisfy 
it is the analytic continuation of the phase space through the threshold: 



2A:^(m) 



above threshold 



p.^.(m) = <( ^ " (38) 

m ^ ^i' below threshold 

It is the basic idea of the so called "Platte model" which has been used to 
analyse the ao(980)-decay into r^vr and KK states, the resonance being very 
close to the KK decay threshold. Similarly, this effect is also present in the 
ii'i(1270)-decays into Kp and K*tt channels with the resonance being at the 
threshold of Kp. This is not so relevant for the iri(1400)-decays where the 
resonance is far above the thresholds. 

The effect of the isobar width. The peculiarity of the i^i(1270) with respect to 
ao(980) case is that the two-body final state includes one unstable particle, the 
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isobar V {V = p or K*). In order to take into account the width of the isobar, 
it is logical to integrate the three-body phase space over the Breit-Wigner of 
the isobar: 



r /27r 

ki{m) = I ki{m,mv) i^^rfmy (39) 

(My - my)2 + '-f 



where ki{m,mv) has its non-relativistic expression 



4 

7 



kAm^my) = \ [m — my—nip) (40) 

V my + mp 

The infinite upper limit in Eq. (39) corresponds to the analytical continuation 
of ki below the threshold for my > m — mp. 

As an approximation to this integral, Nauenberg and Pais proposed to use the 
complex mass of the isobar. My — )• My — iTy /2, in the expression of the momentum 
ki{m,my). These two prescriptions lead to a complex phase space, defined as 



, , 2ki(m) ^ 2 2Mymp f ,^ Vy , ^ 

Piji^) = _ ^ij = —\ I ^ [m- My -mp + i— ] 6ij (41) 



m m]l My + mp \ 2 

where P {P = K or tt) is the final state pseudoscalar meson in the quasi-two-body 
decay. According to us this prescription of using a complex mass is not satisfactory 
for the p and K*, especially for fCi(1270) — Kp. Indeed, we found by direct 
integration of Eq. (39) that the results are quite different from the ones obtained 
using Eq. (41), especially the real part of Pijim) which corresponds to the real 
phase space in the ^^1(1270) — > Kp case (see Fig. 3). The same observation was 
formulated by Frazer and Hendry [37] when the paper of Nauenberg and Pais was 
published. They pointed out that this approximation is valid only for the very 
narrow resonances. The failure of this approach is very worrying since it is basic 
for the whole analysis of Daum et al.. In order to cure this problem, we formulate 
the following assumption: as explained below, instead of identifying the i^'-matrix 
couplings themselves we assume that it is the product of the couplings squared 
and the phase space which is given in a correct way by the experiment, at least 
approximately. 



^For the relativistic phase space Eq. (21) no longer defines a real X-matrix in the physical region. 
The reason is that the relativistic momentum does not remain imaginary below the threshold due to 
an additional complex branch point oc yjim? — (my — mp)'^. Therefore Nauenberg and Pais in Ref. [36] 
restricted to non-relativistic case. 
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m [GeV/c^] m [GeV/c^] 



Figure 3: Dependence of the phase space factor on the mass of the decaying resonance 
m for the K*7t (top) and Kp (bottom) channels. For comparison, pij is calculated using 
the proper analytic continuation, Eq. (39), (blue) and the approximation of Nauenberg 
and Pais, Eq. (41), (red). The difference between two approaches for Re{p22) turns out 
to be significant for the Kp channel. 



• The problem of the P- and D-waves 

In addition, the prescription of Nauenberg and Pais has not been established for the 
P- and D-waves. We do not know what has been done exactly by Daum et al. to 
treat these waves. On the other hand, such waves are to be included in the analysis, 
especially the ktt in the P-wave is very important. Since we are not able to redo 
the analysis by Daum et al. we use the couplings to ii'o(1430)7r channel refitted by 
BABAR collaboration [16]. They include a centrifugal barrier factor depending on 
the complex momentum which is defined by Eq. (41) ^ . However, there is a new 
following problem here. The approximation of BABAR for the centrifugal barrier 
factor is not an approximation to the integral 



1 + kf{m,mv)P? 



(Mv 



(42) 



^According to private communications. 
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which gives a positive real part while the approximation gives a negative one. This 
contradiction can be masked by the normalization of the centrifugal barrier factor 
at the peak. However, this is obviously not a satisfactory solution. 

• The diagonalization of the mass matrix and corresponding rotation of 
the K-matrix couplings into physical couplings 

In several cases we have to deal not with the i^'-matrix couplings but with Breit- 
Wigner parametrization of the intermediate resonances. This is the case, for exam- 
ple, in our calculation of the jT-function. This is also the case of the Dalitz plot 
analyses such as the one of the Belle collaboration [15]. Then the relevant couplings 
are slightly different from those of the i^-matrix. As stated before, they are ob- 
tained from the latter by a complex rotation. Indeed, to pass to the physical states 
we have to diagonalize the mass matrix of the states in the i^-matrix formalism. 
This diagonalization can be performed by a complex orthogonal matrix. This rota- 
tion is complex because of the non-diagonal elements of the imaginary part of the 
mass matrix. The complex rotation angle (which depends on the energy) has both 
real and imaginary parts which are found to be of the order 10° (this result was 
obtained by explicit diagonalization of the mass matrix). As a consequence, this 
rotation affects the couplings: the rotation makes the couplings of the Breit-Wigner 
somewhat different from the ones of the real i^-matrix. The magnitudes of the new 
couplings are different and phases appear. We found that the largest couplings (i.e. 
considering the dominant decay channels, ii'i(1270) — )■ Kp and 7^1(1400) — K*ti) 
are slightly affected and acquire small phases. On the other hand for the smallest 
couplings (i^'i(1270) — > K*tt and 7^1(1400) — )■ Kp) the rotation effects are more 
important. In practical calculations of for the present moment we have neglected 
these effects so that we use directly the couplings obtained from the ^Pq mode ^. 

• Relative signs and "offset" phases. 

It appears that the phases of the amplitudes, deduced from the experimental K- 
matrix analysis are not exactly what is observed: this is a phenomenon of so-called 
"offset" phases. The Kp channel was found to have an additional unexplained phase 
of 30° [8] relative to the {K*tt)s which was set as a reference one. For the ktt channel 
the discrepancy reaches 90°. 

Another problem is that we are not able to establish the complete relation between 
the phase conventions of Daum et al. and quark model ones since the paper of 
ACCMOR collaboration is not detailed enough. 

^For more details, see the Appendix A 
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5 Numerical results 



Let us summarize our final prescriptions we use for tlie calculation of the partial widths 
and for the further extraction of our theoretical model parameters from the experimental 
measurements. Our basic approach is to use partial widths at the peak on both, theoretical 
and experimental, sides. We abandon the idea of using the branching fractions and the 
total i^i-widths for the comparison with our predictions. 

1. For the theoretical prediction, in order to take into account the isobar width effects 
in our theoretical prediction of the partial widths T^^*"*^, the amplitudes (9) squared 
are integrated over the invariant mass of the isobar: 



r^ii =8^ / —T7 \Ai{Ki^VP)\ ^dmv (43) 

Jm™" ^^-^Ki (My - my)2 + '-f 

Note that since we consider the widths at the peak there is no integration over 
the Ki invariant mass unlike what is done in several theoretical papers (e.g. see 
Ref. [28]). Moreover, one can notice that the integration over the mass of the isobar 
is one within the correct physical region restricted by the corresponding physical 
bound of the two-body decay (i.e. we use the real phase space). 

2. For the experimental input, we make the simple assumption that the partial widths, 
calculated from the i^'-matrix couplings at the peak according to Eq. (44), are 
correct, although the complex phase space a la Nauenberg and Pais (41) might be 
not correct (i.e. what we measure by fitting data, is always the combination like 
fa'(b')i ^ Piji^) which are assumed to be extracted correctly). Therefore, we use the 
/T-matrix couplings and the real part of the complex phase space a la Nauenberg 
and Pais in order to extract the experimental values of the partial widths 

r^xt = ^fa'm^Re[p^AM,eak)] (44) 

3. We calculate this partial width according to Eq. (24) also for the P (L = 1) and D- 
waves (L = 2), assuming that the i^T-matrix couplings / contain the barrier factors 
B[{m) that are properly normalized at the peak: 



/a(6)i("T')|p r)_wavcs ~ •^'»(*)* uL 



BHm) 
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. n L/2 (45) 



1 + kf{m)R^ 



where = 25 GeV~^ [16]. This assumption seems to be correct since it leads to 
the calculated branching ratios that are very close to the ones announced in the 
paper by Daum et ai. In any case we avoid as much as possible to rely on the 
experimental data on ifi(1270) — )■ Kp and the Z)-wave of Ki(1270) — t- K*tt and we 
trust our theoretical prediction. 



5.1 Fit of parameters 7 and Ok^ 

In order to extract our phenomenological parameters, the quark-pair-creation constant 7 
and Ki mixing angle, we do a fit using the method of least squares. As an experimental in- 
put we use the partial widths (namely, r^"f from Table 4) only of the following processes: 
/s:i(1270) ^ {K*7i)s, Ki{UOO) {K*7f)s, Ki{UOO) {Kp)s, which are assumed to be 
Gaussian distributed with mean r^^'"^^(7, 9xi) and known variance a^peak. The Z)- waves 

are not taken into account in our fit. Moreover, the dominant channel ii'i(1270) — )■ Kp 
due to the dangerous threshold and phase space effects is avoided since the narrow width 
approximation can be incorrect for the decays near the threshold and here the width 
effects can play a significant role. 

Then, the likelihood function is constructed as a sum of squares 

3 /-ppeak _ -pQPCM / n \\2 

X\l.e,,) = -21nL(7,^Kj = E " 2 ^ (46) 

i=i '^rP""'= 

In order to find the unknown parameter O^i the function is minimized, or equiva- 
lently the likelihood function L{6ki) maximized. The minimization of the gives the 
minimal value Xmm = 0-61 and the estimators 7 = 4.0 and 9ki = 59°. 

The covariance matrix for the estimators Vij = cov[^i,C,j] can be found from 



1 d'x' 



(47) 



Thus one obtains 



<yl \ / 0.29 0.99 



-"-"-'^ic,;. <-j-(o:99 107,0 j (^«) 

where the diagonal elements give the variances a? and a? . Finally, one finds the fitted 



values of the quark-pair-creation constant and Ki mixing angle: 

7 ~ 4.0 ± 0.5, Ok^ ~ (59 ± 10)° (49) 

Taking for granted that our theory is correct, one is now interested in the quality of 
the agreement between data and various realizations of the theory, determined by the 
set of parameters, namely {-^.Ok^}. For metrological purposes one should attempt to 
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estimate as best as possible the complete set of parameters {7,6'^-^}. In this case we use 
the offset-corrected [38]: 

Ax\l,eKj = x\l,OK,)-xLn (50) 
where Xmm absolute minimum value of the function of Eq. (46) which is obtained 

when letting our model parameters free to vary. The minimum value of Ax^ is zero by 
construction. Here one has to notice, that this absolute minimum does not correspond 
to a unique choice of the model parameters. This is due to the fact that the theoretical 
predictions used in the analysis are affected by important theoretical systematical errors. 
Since these systematics are restricted in the allowed regions there is always a multi- 
dimensional degeneracy for any value of x^- However, since in our analysis there are 
only two model parameters, our predictions for {7,6*^-1} are not affected by any other 
theoretical predictions. 

A necessary condition is that the confidence level (CL) constructed from Ax^(7,^i^i) 
provides correct coverage is that the CL interval for {7,6'/^^} covers the true param- 
eter value with a frequency 1-CL if the measurements were repeated many times. The 
corresponding CL intervals for the confidence level of CL=68% are shown in Fig. 4. 

5.2 Model predictions for partial widths 

Now, we can make systematic predictions for various processes. First, it is very useful to 
check our result for the quark-pair-creation constant 7 prediction with the much better 
studied 61 — )■ {ujn)s and bi — > (wtt)/) decays which depend only on 7. One can see 
from Fig. 5 that our estimation for 7, determined from the i^'i-decays (49), is in a good 
agreement with the one extracted from the bi — )■ utt decay. Moreover, the extracted D/S 
ratio of the partial amplitudes is very well predicted and coincides with the measured 
value including the sign: 

{Ad/As)qpcm ^ 0.28 (51) 

while the experiment [39] gives: 

(An/AsUp = 0.277 ± 0.027 (52) 

Note that the Belle collaboration omits the waves in the B — ?■ J/ipKi analysis. This 
could be of consequence, since the Dalitz plot should be appreciably different according 
to our calculation (see our discussion in the end of subsubsection 5.3.3) 

^^In statistics, a confidence level interval is a particular kind of interval estimate of a fitted parameter 
and is used to indicate the reliability of an estimate. It is an observed interval (i.e. it is calculated from 
the observations) , in principle different from sample to sample, that frequently includes the parameter of 
interest, if the experiment is repeated. How frequently the observed interval contains the parameter is 
determined by the confidence level. 

^^One has to point out that the branching ratio of hi — > ujtt has not been measured precisely. However, 
the WTT is considered to be the dominant decay mode [39], so that we assume B{bi — >■ wtt) ~ 100%. 
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Figure 4: distributions for the fitted parameters, Ki mixing angle 9ki and quark- 
pair-creation constant 7 (left), with the confidence level intervals that determine how 
frequently the observed interval contains the parameters (right). 

To summarize, we give in Table 5 our predictions for the S'-wave partial widths of the 
strong interaction decays of the .fri-mesons, using the fitted values of 7 and 9ki ■ One can 
see that the agreement is satisfactory except for the ii'i(1270) — )■ Kp channel. This is 
not unexpected in view of the particular difficulties of the experimental treatment in this 
decay as explained in the previous section (recall especially that the drawback of using 
the phase space formula of Nauenberg and Pais is crucial in this case) . 



Decay channel i 


MeV/c2 


Mev/c^ 


iri(1270) -> {K*7r)s 


31 


28±26 


i^i(1270) ^ {Kp)s 


61 


122±28 


^^1(1400) ^ {K*Ti)s 


209 


211±59 


K,{im) ^ {Kp)s 


1 


20±25 



Table 5: Theoretical predictions for the partial decay widths, calculated using the fitted 
parameters 7 = 4.0 and 9ki = 59° and compared to the experimental partial values of 
widths at the peak (see Table 4). 

As for the D-waves in the i^Ti-decays, our impression is that they are poorly determined 
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Figure 5: QPCM constraints for the quark-pair-creation constant 7 and the Ki mixing 
angle Ok^ obtained from the fitted partial decay widths at the peak, calculated using 
the i^-matrix couplings (Table 4). The cross indicates the optimal values of 7 and 9ki 
extracted from the fit. 



experimentally. Our prediction (r(_ft'i(1270) -> {K*-k)d) — 3 MeV/c^) lies below the 
experimental numbers: the couplings for the D-waves are not given in the paper by Daum 
et ai. Tentatively they were re-fitted by the BABAR collaboration [16] from which we 
deduce the partial width r(/s:i(1270) ^ {K*7t)d) = (34 ± 3) MeV/c^. Here one has 
to notice that the errors of the re-fitted parameters are surprisingly small, as the ones 
obtained by Daum et ai. 

5.3 Prediction of signs of decay amplitudes and the "offset" 
phase issue 

Let us recall that, at least for the determination of the photon polarization parameter 
as described in our paper [3], our goal is to calculate the jT-function (1) which describes 
the full three-body Ki — )■ Kttt^ decay. As explained, we need in fact the expression 
Im[n- {J' X J'*)] which depends crucially on the relative phases of the V — )■ PP couplings 
and the Ki — )■ VP form factors (see Eqs. (22)-(27) in Ref. [3] for the definition). These 
quantities are directly related to the two-body decay amplitudes, calculated by using the 
quark model. The phases of these amplitudes do not make sense by themselves but only 
in the product of two amplitudes of the subsequent processes which describe the final 
three-body decay Ki — )■ Ktth. Then, the relative signs are observable quantitities, that 
can also be determined from any careful experimental study of.the Ki decays. We define 
the relative phases for two Ki — )■ Kirn amplitudes of various partial waves via different 
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intermediate isobar states (i.e. {K*7!-)s, {K*tx)d, {Kp)s). Standardly, the reference partial 
wave is chosen to be the 5*- wave of i^*7r. For instance, the relative phase of the Ki — )■ 
Kp — )■ i^TTTTT channel is defined as: 

^ r AsjK, ^ Kp) xApjp^ 7i7c) 
" ~ ""^^ [As{K^ K*n) X Ap{K* Kn) 

One has to note, that the total relative phase which is contained in the j7-function contains 
of course complex the phase of the denominator of Breit-Wigner of the isobar. For the 
conventions necessary to define 6p we refer to Appendix. 

6p is independent of the conventional phase factors of the meson states (e.g. meson 
wave functions ^^). In the ^Pq model each decay amplitude is real with suitable conventions 
of the wave functions and by factorization of spherical harmonics. Then in the quark 
model 6p is real. This is due to specific properties of the transition operator. 



(53) 



5.3.1 Sign of the D/S ratio 



The simplest prediction is the one concerning the D/S ratio in the bi — t- uir and ai — )■ pvr 
decays. Indeed, this sign depends only on the well known standard conventions. It is then 
striking that all the signs are correctly predicted by the model. In the case of bi and Oi 
these signs are well measured and given in PDG. For the Ki — )• K*7t channel the signs 
are not given by Daum et al. in [8]. However, we can read the relative phase for i^'i(1270) 
from Fig. (13) in Ref. [8] which is positive {fb^/ fbi > 0), while for fri(1400) we have to 
rely on the analysis of BABAR because it is not possible to fix it from the figure since the 
D-wave is too weak overwhelmed compared to the D-wave of 7^^1(1270) {fa^/fai < 0). 

In the paper of Gronau et al. [1,2] the D/S" phase for K*tt is given as 6d/s = (260±20)°. 
We believe that the authors were misled by incorrect interpretation of Fig. (13) (bottom- 
right) in [8]: the plotted phase indeed peaks at 260° at Mk-k-k ~ 1-4 GeV/c^ . But this 
is not the phase we are looking for since it contains the phase from the Breit-Wigner of 
ii'i(1270) which is dominating over the i^i(1400) contribution and gives an additional 
phase of approximately 90°. Hence, the phase we are interested in must be read as 
^D/s ~ (260 — 90)° ~ 180°. We must stress the following subtle point: the plotted 
phase is the difference of the phases of the D-wave strongly dominated by i('i(1270) 
and the one of the S'-wave which includes large contributions of both resonances. As a 
consequence, paradoxically, there appears a bump in the D-wave phase diagram, peaked 
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In the QPCM, 5p can be calculated from 



Sp cx arg 



arg 



^(Ki)^(K)*^(-n)*^(-n)* 



1 



(54) 



what impUcs that the relative phase of the total amplitudes is real (i.e. = or tt) and does not depend 
on the separate complex phases of the meson wave functions. 
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at Mk-k-k ~ (1-3 — 1.4) GeV/c^ which is essentially determined by the tail of the Breit- 
Wigner of i^i(1270). We checked this conclusion by explicit calculation of the amphtudes 
using the iiT-matrix couplings (see Fig. 6). 



240 

230 

M 220 
a 
-a 

- 210 
200 
190 
180 

1.0 1.1 1.2 1.3 1.4 1.5 1.6 

m [GeV/c^] 

Figure 6: The D-wave phase relative to the S'-wave of K*tt, calculated using the iC- matrix 
couplings. One can see a bump at Mk-k-k ~ (1-3 — 1.4) GeV/c^. 

5.3.2 Relative sign of the Kp/K*Tx couplings 

We study the real phase (i.e. the relative sign) of the i^'i(1270) K*!^ and i^ri(1270) -)• 
Kp amplitudes, which plays important role in the A-y determination using the cu-method 
(due to the strong dependence on the phase of the interference term Im[n ■ [J x J*)])- 
Indeed, the odd moments of u change their sign if one changes the relative sign between 
the -> K^p^ — )■ A"+7r"7r+ and — )■ _ft'°*7r+ — )■ A"+7r"7r+ amplitudes. One has to 
notice that in this case this phase can be hardly extracted from the i^-matrix analysis 
by Daum et al. due to some unknown conventions (in particular, the order of particles 
what is significant for the determination of the couplings signs). We then rely on the 
recent analysis by the Belle collaboration of the B — )■ J /ip{ip')K'K'K decay which gives 
more explicit explanation of the conventions. 

Here we summarise what is new in the Belle B — )■ J /ip{ip')K'K'K paper [15]. First we 
will list up the general conclusions of this paper and then, discuss some details of the 
Dalitz plot shown in this paper, which provides important information to our work. 

5.3.3 General conclusions of the study of — t- J/^fCvrTr by the Belle collabo- 
ration 

This paper, in principle, focuses on the measurement of the branching ratios of B^ — )■ 
J/ipK^Ti^Ti^ and B^ — tp'K'^Ti^n^ . Since the Ktctc final state comes from various 
resonances, -K'res; this analysis provides information of the Kj-^s — > Kittt strong decays. 
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Since the K^^s = -ft'i(1270) turned out to be a prominent component (for both J/ip and 
ip'), some detailed study of ^1(1270) — )• Kim has been done: 

• The Dahtz plot for the three-body decays is shown. We discuss more details on this 
later. 

• The intermediated two-body decay branching ratios have been re-determined (see 
Table 6). The branching ratios for the dominant decay modes, ii'i(1270) — )■ Kp and 
ii'i(1270) — > K*TT, are found to be slightly different from the previous measurements 
(PDG), although they are still in accordance within several standard deviations. On 
the other hand, the ii'i(1270) — )■ -K'Q(1430)7r channel, which was supposed to have a 
large branching fraction (-B(i^i(1270) ^ K*{U30)7r) = (28 ±4)%) according to the 
previous measurements [8,39], was found to have a significantly smaller contribution 
of the order of 2% (see Table 6). 

• In addition, by floating the mass and width of the i^'i(1270) in an additional fit of the 
5+ J/tpK+TT+TT- data, a smaller mass of (1248.1 ± 3.3(stat) ± 1.4(syst)) MeV/c^ 
and larger width (119.5 ± 5.2(stat) it 6.7(syst)) MeV/c^ were measured for the 
i('i(1270). Of course, there is a correlation between the fact that the "scalar -|-7r" 
component becomes much smaller and the fact that the K*'n' and Kp contributions 
become larger (see Table 6). 

Here we want to draw attention of the reader to the conceptual difficulties raised by 
the definition of the ii'i(1270)-width. In the Fit 1 the Ki width is the one given by 
PDG while in the Fit 2 the width was treated as a free parameter. Due to the threshold 
effect one should not expect that the width measured by the Belle collaboration from 
the Breit-Wigner denominator at the peak should coincide with the one defined by PDG, 
although it should be much larger. One observes that the floated width is larger than the 
PDG value but it is still much smaller than 200 MeV/c^ as we would expect from the 
calculation using the i^'- matrix formalism (see Table 3). 

One has to point out that the D-waves are not taken into account in the master formula 
of Belle. On the other hand, we found from the theoretical study that the D-waye of K*7i 
can have a small but non-negligible effect. In principle, there are two bumps due the 
presence of the D-wave, but it is found that the one located in the intersection region of 
the Mxn ~ Mx* and M^rTr ~ Mp on the Dalitz plot is masked by the dominating peak of 
p. Using a Monte-Carlo simulation, we observed a second small but non-negligible bump 
at low Mtt^ (see Fig. 7 in the center). 

5.3.4 Dalitz analysis 

In [15], the Dalitz plots for iri(1270) — t- K-kh is shown in the three variable planes, 
M2(K+7r+7r-), M\K+n') and M^{n+TT-). On the Dalitz plot in the M\Kn)-M\TnT) 
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Decay mode 


PDG (%) 


Fit 1 (%) 


Fit 2 (%) 


Kp 


42 ±6 


57.3 ±3.5 


58.4 ±4.3 


K*TT 


16 ±5 


26.0 ±2.1 


17.1 ±2.3 




28 ±4 


1.90 ±0.66 


2.01 ±0.64 



Table 6: The fitted branching ratios of the i^'i-decays measured by the Belle collaboration 
in the analysis of 5 — )■ J/ipKirn decay [15]. 

plane, a strong interference effect between Ki — )■ K*7r and i^ri(1270) — )■ Kp is observed 
(see Fig. 7). In particular, it is pointed out that the weakening of the Kp in the region 
of M{K7t) > Mk*{s92) is originated from the interference of the Kp and K*tt amplitudes. 
Here we will attempt to study the real phase (in another word, the relative sign) of the 
Ki — )• if*7r and -ft'i(1270) — )■ Kp amplitudes using this Dalitz plot, to check our theoretical 
prediction. Indeed, as we will see later-on, in a forthcoming paper, this information of 
the phase has an important consequence on our determination. 




Figure 7: Dalitz plots of — )■ i^'j*"(1270)7 — > K^ti^'k^'j, measured by the Belle col- 
laboration [15] (left) and MC simulated for the "offset" phase equal to (center) and vr 
(right) of the Kp channel relative to {K*'k)s. The "correct" phase Sp = corresponds to 
our quark model prediction. 

5.3.5 Determining the relative sign of the Kp/K*^ amplitudes 

In this section, we demonstrate how the relative phase between the Kp/K*7r amplitudes 
can be determined from the Dalitz plot. 

In [15], the full amplitude of Ki three-body decays is defined as 

\M{sKi,SK*,Sp)\^ = \aK'AK*{sK^,SK*) + apAp{sKi,Sp)f (55) 
where the coefficients aK*,p represent the strong decay of Ki — > Kttit through K*, p 
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intermediate states. The amplitudes Ak*,p are defined as 



where Pi is the breakup momentum of Pi or V in the i^'i reference frame and is the 
angle between the momenta of Pi and Pk in the V reference frame, which can be expressed 
in terms of Sk,, Sij, Sik 

Compared to the obtained Dalitz plot, we can determine the coefficients aK*,p in- 
cluding the relative phase between them. The obtained result by the Belle collaboration 
yields [15]: 



\aK* I = 0.962 ± 0.058 ± 0.176, |ap| = 1.813 ± 0.090 ± 0.243 

(5p = arg(ap/axO = -(43.8 ± 4.0 ± 7.3)° (57) 

Formula (55) can be written in the following general form factorizing out the phase: 



\M{sK,,SK*,Sp)\'^ = Co{sK,,SK',Sp) + Ci{sk,, Sr* , S p) COS Sp + C2{sk,, Sk* , Sp) sm 5p (58) 

where Cj(s, sk-k, Stttt) are the known functions, expressed in terms of various combinations 
of the real and imaginary parts of \aK*\AK*{sKi, sk*) and \ap\Ap{sK-,, Sp). So, in order to 
establish the correspondence between our parametrization of |A^p (|i/P in our case) one 
can compare the relative signs of the cos5p and sin5p coefficients, Ci 2, on the Dalitz plot. 
Direct numerical calculation shows that 

sign {cr""') = sign (cf ^"^) , sign (cr"^') = -sign (cf ^"^) (59) 

5.3.6 The issues of complex "offset" phases 

In principle, the QPCM predicts real Ki — )■ VP amplitudes, without any complex phases. 
This should correspond to the i^-matrix couplings. The complex rotation of the iiT-matrix 
states to the physical states should introduce complex phases but we found by explicit 
calculation that the imaginary part of the rotation angle is small: 

iPa'^aP>^ - 10° (60) 

However, the Belle collaboration measured a sizebly larger imaginary relative phase 
(i.e. Eq. (57)) of 6p ~ —44°. We recall also that Daum et al. measured a non-zero phase 

^•^One has to notice that the D-wave amphtude is not taken into account in this parametrization and 
that the last factor in Eq. (56) corresponds to the ^-wave. 
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of the order of 30°. Similar value was found in the reanalysis of the ACCMOR data by 
the BABAR collaboration: 6p = -31° [16]. 

There is no explanation of this complex phase in a definite theoretical model: neither in 
the ^Pq quark model nor in the most general quasi-two-body i^-matrix approach. Indeed, 
the "offset" phase which is introduced in the analysis by Daum et al. depends only on the 
decay channel and is the same for the lower and upper resonances. The general production 
amplitude for each channel in the reaction K~p — )■ {K^-k^-k^) is written as [8, 16] 



where the factor — represents the propagation and the decay of the i^'i-resonance. 

The last factor Pj describes the resonance production which can be in principle complex 
(indeed, one finds in [8] that there is a non-zero relative phase between the production 
couplings of two i^i-resonances). From Eq. (61) it is obvious the "offset" phase 5i can not 
be ascribed to either the resonance decay or production amplitude. 

This puzzling situation must not be ignored and has to be studied more carefully. 
In the present, we use the model prediction for the j7-function as it is with pure real 
couplings. On the other hand, to adopt pragmatic attitude we explore the effect of the 
introducing this additional "offset" phase 5p = — (J^eiie ^^^q calculation of the j7-function 
and the estimation of the theoretical uncertainty of A-y. 

5.4 The issue of the kit channel 

The PDG assigns a large branching ratio B{Ki{1270) K*{U30)7r) = (28 ± 4)%. It is 
extracted as all the branching ratios, from the ACCMOR data and analysis. However, 
this interpretation is dubious. The original ACCMOR measurement shows indeed a clear, 
strongly coupled peak in the "scalar + vr" channel around the mass Mk-k-k ~ 1270 MeV/c^. 
However, it is not at all claimed that the scalar is ii'Q(1430); it is treated as a lower 
and much broader scalar meson (M ~ 1.25 GeV/c^, F ~ 600 MeV/c^); or could be a 
continuum {KT\:)s-wa.wc according to [40]. 

The -ft'o(1430)-meson is the scalar orbitally excited state of kaon which has the mass 
Mk*{iazo) = (1425±50) MeV/c^ and width F^. (1430) = (270±80) MeV/c^ [39]. According 
to quark models, the constituent quarks are in the ^Pq state. In order to estimate the 
i^'Q(1430)7r contribution we use QPCM to calculate the P-wave amplitude for the decays 
7^1(1270) K*{im)7r. One can see from Fig. 8 that Ap{Ki{1270) i\:*(1430)7r) 
is strongly suppressed compared to As{Ki{1270) — )■ K*7i). Moreover, there is also a 
suppression due to the phase space. Finally, after the integration over the phase space for 
y/SK-K within the allowed physical range [niK + rri^; Mi^j(i27o) — 171^,], we predict that 




(61) 



B{Ki{1270) K*{U30)7i) 
B{Ki{1270) K*{892)7t) 



< 0.01% 



(62) 
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in blatant contradiction with the PDG entry. 

What is most striking is that indeed, the Belle collaboration finds B{Ki{1270) — )■ 
ii'g (1430)7r) ^ 2% (see Table 6); it is very small as we predict. They did not find 
any other "lower scalar+vr" component in the i^i-decay: the B missing with respect to 
ACCMOR seems to be filled by an enlargement of Kp. Therefore, in our analysis, we do 
not include the i^i(1270) — )■ i^o(1430)7r channel. Neither do we include any other possible 
scalar in the presented results. However, to take into account the contrary conclusions 
of ACCMOR, we keep in mind the possibility that there is some significant portion of 
the branching ratio carried by a very wide scalar meson, different from the i^'Q(1430), 
such as the low lying state ii'Q(800) (also called k) [41]. Note that such state is most 
probably not a qq state and therefore the decay into ktt can not be estimated within our 
theoretical model. Such contribution has not been tested explicitly in the analysis by the 
Belle collaboration. 

Let us mention two other relevant facts: on the one hand the non-strange counter 
part of k(800), cr, is found with sizable branching ratio in the decay of ai(1260) in the an 
state. On the other hand it is surprising, as noticed by Daum et ai, that there is no kit 
channel in the i^i(1400)-decay. 

7 
6 
5 

._4 
~ 3 

2 

1 



0.4 0.6 0.8 1.0 1.2 

.Vk, [GeVVcl 

Figure 8: \As{Ki{1270) K*7t)\'^ (red) and \Ap{Ki{1270) -> K*{im)7T)\'^ (blue) for 
ski = ^1-1 (1270)- "^^^ -^1 i^ixing angle Ok^ is taken to be 60°. 

6 Conclusions 

Let us now summarise the main conclusions of the present work, and sketch some prospects 
for progress regarding theory as well as experience. 

Not only the strong decay pattern of ii'i-mesons is quite complex, but, not surpris- 
ingly, it is then difficult to analyze the whole system experimentally. In spite of many 
efforts, we have found that much information is lacking, and that certain weaknesses 
may be suspected in various analyses. In lack of more fundamental treatments, we have 
recoursed to the quark model approach to explain and complement the experimental re- 
sults. The quark model, although approximate, is the basis of our whole understanding 
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of spectroscopy. The ^Pq model for decays presents the advantage of handhng in a simple 
way the whole set of L = 1 decays On the other hand, experimental input is still 
required to fix necessary parameters, for instance the Ki mixing angle. 

Our predictions for the specific strong decays under concern, i.e. L = 1 states de- 
caying to VP states, can be evaluated by comparing to data where available. On the 
whole, our conclusion is very encouraging. In addition to the known fact that a certain 
mixing is able to explain the pattern of VP decays, the model explains detailed features 
which are quite outside an SU{3) symmetry approach, and require a dynamical approach. 
This is the case, for instance, the D/F ratio of octet couplings and the D/S ratios in 
magnitude and phase; it is an achievement of the model that all the observables phases 
are correctly predicted. Another typically dynamical prediction is that the decays to qq 
scalar+pseudoscalar should be very small. It is in agreement with the recent observations 
by Belle [15], but it does not exclude a large contribution of non qq scalars, which could 
then explain the observations of ACCMOR. In any case, it strengthens the conclusion that 
the K channel observed by ACCMOR is not the K^{1430) (presumed qq) as tabulated in 
the PDG tables. 

We beheve that the Ki system deserves further investigation because it has revealed 
interesting in various aspects. Indeed, it also presents unexplained features in the standard 
domain of spectroscopy, i.e. the mixing angle and the mass splitting. The mixing of the 
two states offers the possibility to explain the remarkable pattern of K*n/Kp decays, but 
the angle is not a theoretical prediction. In fact, in a potential model, spin-orbit forces 
generate a mixing, but it is not the one which is observed. As explained in the text, loop 
effects would also generate a mixing effect, but it cannot be calculated. It must also be 
noted that the mixing does not explain why the ktt channel if present in i^'i(1270) is absent 
in i^'i(1400): quite on the contrary, as explained in the paper of Daum et al, one would 
expect the mixing to generate a coupling from the Ka component. Apart from mixing, 
the predictions of potential model also fail to explain the splitting of the two states - it 
is predicted much too small by the model of Godfrey and Isgur, which can be estimated 
to be the most trustable. These facts show that our knowledge of spectroscopy is not yet 
satisfactory even for apparently well identified, low-lying qq states. Understanding these 
facts then justify further studies. 

We try to say something about possible improvements: 
1 Improvement of the theoretical treatment 

It is important to recall that there is no fundamental theoretical treatment of such prob- 
lems and that quark models, on which our theoretical model is based, although much 
valuable, contain essential approximations, i.e. ones that cannot be improved systemati- 
cally. This holds in two respects: potential models are of course essentially approximate, 

"'^^In fact it is supported by a much larger set of experimental tests. 
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even with relativistic improvements as included in the model of Godfrey and Isgur, but this 
is also true of the quark-pair-creation decay model itself, which, presently, is essentially 
non-relativistic. The center of mass motion of the hadrons is not treated relativistically. 
Progress is desirable in this direction. 

2 Prospects of improvement of experimental knowledge 

At present, further progress could come mainly from a better and more complete de- 
termination of the magnitudes and the phases of the various couplings by experiments. 
Certainly, the old experiments with production of Ki by strong interaction scattering, 
as the ones of SLAC and ACCMOR, have much larger statistics for decays involving Ki 
than present B factories . Yet there is little prospect of them being redone, and they 
have also their own weakness in the fact that the production process is complex. On 
the other hand, there is the hope that new detailed studies could be made in B and 
r-decays. Encouraging examples have been coming from both BABAR and Belle, such 
as r — )■ KiUr [42] and B — )■ Kiip [15]. In fact, a distribution in an additional angle may 
also help to improve the analysis. For example, a new study of B ^ Kiip with angular 
analysis could yield directly the crucial quantity Im{n- {J x J*)) up to a multiplicative 
constant [3]. The analyses could be guided by our semi-theoretical and approximate in- 
vestigation, which, for instance, emphasizes the need to take into account D waves, not 
included in the present Belle analysis of i? — ?■ Kiip. 
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A Re-interpreting the ACCMOR result in terms of 
"physical states" 

In order to determine our model parameters and the Ki mixing angle from comparison of 
the predicted partial decay widths of the i^i-meson decays into the dominant i^*7r and Kp 
channels with the measured experimental values, we use the fitted iiT-matrix parameters 
extracted by Daum et al. from Ref. [8] (see Table 8). 

Using the definition of the Ki mixing by Daum et al. (which is different from (2) that 
we use by signs): 



|iri(1400)) = |i^i^)cos^;^, + |iriB)sin0^, 
|iri(1270)) = -\K^a) ^i^Ok, + \K^b) cosOk, 



(63) 



the dominant S'-wave i^-matrix couplings of the Kis to the states K*7t (channel 1) and 
Kp (channel 2) are given as [8] 



fa'l = 


^7+ COS 6k^ + \ 


fb'l = 


-^7+sin^Xi + 


fa' 2 = 


^7+ cos 9k, - \ 


fb'2 = 


-^7+sin^/<i - 





7_ sin 



7_ cos 

(64) 

7_ sin 9ki 



7„ cos 6k, 



where 7+ and 7_ are the reduced SU{3) couplings for Kia (F-type) and Kib (D-type) 
respectively. Their fitted experimental values are given in Table 7. The indices a' and b' 
denote the upper and lower Ki resonances. 





rua' , GeV/c^ 


mf,r , GeV/c^ 


7+ 


7- 


6ki 


1.4±0.02 


1.17±0.02 


0.78±0.1 


0.54±0.1 


64°±8° 



Table 7: Fitted i^'-matrix pole masses, S'-wave reduced SU (3) couplings and mixing angle 
for KiA (F-type) and Kib (-D-type), taken from Ref. [8] (low t data). The indices a' and 
b' denote the upper and lower Ki resonances. 

Using the experimental values of the i^-matrix couplings from Table 8 and performing 
the diagonalization of the complex mass matrix (31), we observed that 

• The variation of the absolute values and phases of the new rotated physical couplings 
{faphi, fbphi} around the masses at the peak of Breit-Wigner (i.e. m ~1.27 GeV/c^ 
and 1.4 GeV/c^) turn out to be small (see Fig. 9). 
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fa'l 


fb'l 


fa' 2 


fb'2 


fa' 3 


fb'3 


0.50±0.07 


-0.19±0.09 


-0.15±0.10 


-0.51±0.06 





0.32 



Table 8: i^-matrix couplings, calculated from Eq. (64) using the fitted parameters from 
Table 7. The indices a' and b' denote the upper and lower Ki resonances decaying into 
K*7i (channel 1) and Kp (channel 2) hadronic states respectively. The coupling to the 
Kq{14:30)71 channel, where i^'o(1430) resonance is supposed to have the mass 1.25 GeV/c^ 
and width 600 MeV/c^, /^s is taken from Ref. [16]. 

• Contribution of the complex phase space for energy below the decay threshold (which 
implies Pij{m) — > i\pij{m)\) is very small for diagonalized physical mass of -ft'i(1400) 
(see Fig. 10). But one observes a threshold effect for 7^1(1270) near m ~ 1.2 GeV/c^. 
However, the mass variation of Mxi{m) around the peak of Breit-Wigner can be 
considered not so significant. 

• One can see from Fig. 10 that, contrary to M/^^(m) dependence, the width r^^(m) 
is a rapidly varying function of the energy m. 

• Non-diagonal elements of the mass matrix (31) are sufficiently small compared to the 
diagonal ones. One can see from Fig. 10 that the difference between the properly 
diagonalized masses and widths (blue/red curves), which are calculated in terms 
of the rotated physical couplings, and the real and imaginary parts of the diago- 
nal elements of (31) (green/orange curves) is insignificant. As a consequence, our 
assumption for the partial widths 

r„p.,(Mpeafc) ^ r,,,(Mpeafc) = rJf^^(Mp,,,) (65) 

seems to be reasonable. This means that we can use the experimental measured K- 
matrix couplings in order to calculate the partial decay widths and fit our moedel 
parameters, namely quark-pair-creation constant 7 and the mixing angle 6ki , which 
can further be used for the J7 function computation. 
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Figure 9: Energy dependence of the physical couphngs (blue). The red lines represent the 
values of the real couplings for the i^'-matrix states, fitted by ACCMOR collaboration [8]. 
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0.30 




Figure 10: Energy dependence of the mass (left) and width (right) of -ft'i(1270) (red, 
orange) and -ft'i(1400) (blue, green). Red and blue curves correspond to the masses and 
total widths of the physical eigenstates, i.e. diagonal mass matrix elements (??) which are 
calculated in terms of the rotated physical couplings. Orange and green curves represent 
the leading diagonal elements of the complex mass matrix (M' — iV /2)^,y (??) in the 
i^'-matrix eigenstate basis. The D-wave contribution is not taken into account due to the 
absence of knowledge of the corresponding couplings. 

B QPCM 

B.l Spacial integrals in QPCM 

For the decay A B + C (see Fig. 1) the spacial integrals are given by 



I^^BC) ^ j d%d%d%d%5{h + k2- kA)5{k2 + h- kB)S{h + h - kc)5{h + k^) 
= ^ j d^ky^{kB-k)ij^^\kB + k)ij^''\-k)i;^^\k) 

(66) 

where ip^s are the normalized Fourier transforms of harmonic oscillator meson wave func- 
tions. The wave functions for the ground (L = 0) and orbitally excited (L = 1) meson 
states are defined as 

yr{k) = \k\Yr$) = {eJ)^ 
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Here Ri is the meson wave function radius and Sm are the A-polarization vectors, defined 
as 

\ /I 




£±1 = T 



1 

71 




(68) 



Performing the integration over k one obtains for the orbitally excited axial-vector 
meson decay into pseudoscalar and vector mesons in the ^-meson reference frame: 



^ ~ 7r5/4 (R\ + RI + Rl,)5/2 [^^"^ • ^ • '^sA^-m • Kb) 



X 



{2R 



11 Q 



RhKR 



A{R\ + RI + Rl) 



exp 



f^B 



R\{R 



B 



8{R 



i?2 



Rl). 



(69) 



Setting ks along z-axis, the integrals become 



AABC) _ 4\/3 R^/^{RbR, 



c. 



1 3/2 



7r5/4 (R 



i?2 



Rb 



B 



A{R 



^B 



R'c) 



X exp 



-k 



Ra{Rb + Rc) 

8(R\ + R% + Rc) 



(ABC) 



AVS RfjRBRc)^'' 
'7rV4 {R\ + Rl + Rlf/-^ 



exp 



-k 



R^jRB + Rc) 

8(R\ + R% + Rc) 



Rc^ 



(70) 



For the vector meson ground state decay into two pseudoscalar mesons the spacial 
integral is 



j{ABC) 



V6 



5/4 



(£m • kc) 



{RaRbRc)'/H^Ra 



Rb + Rc) 



X exp 



{R\ + RI + Rlfl^ 
R\{R 



-k 



Rl) 



C of td2 I p2 I td2\ 



8(i?^ 



(71) 



B.2 Fixing the relative signs for three-body decay 
B.2.1 Clebsch-Gordan coefficients 

As it was emphasized in the text the relative sign of several amplitudes involving various 
intermediate states plays very important role. Therefore the convention of the particle 
order in the Clebsch-Gordan coefficients is very important. For instance, for the case of 
K^7T~TT^ decay that implies that we take the Clebsch-Gordan coefficients defined 
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in the following way: 



{K*\+\Kt) = (1/2, -1/2; 1, l|l/2, 1/2) ' ^ 



3 

2 



{K+ir-lK*") = (1/2, 1/2; -1, 1| - 1/2, 1/2) , ^ ^ 
iK+p'\Kl) = (1/2, 1/2; 1, 0|l/2, 1/2) = 
(vr-7r+|p°) = (l,-l;l,l|l,0): ^ 



(72) 



x/2 

This gives the signs of the amplitudes listed in Table 2. 

B.3 Determination of the relative sign of gK*KTr and Qp^^^^ 

Following the definition in the work of Gronau et ai, the total amplitude of the two 
possible channels is written as 

M^) = ef^^T^yj'>g,.A'\p.^ -P.-Y 

where Typ is the hadronic tensor, parametrized in terms of the form factors /y, hy (or 
equivalently the S and D partial wave amplitudes) 

Now, using the same Clebsch-Gordan coefficients, defined above in Eq. (72), one can 
write the amplitude of the V — )■ PP decay, calculated the general tensor Lorenz-invariant 
form in the vector meson reference frame: 



M{K*^ AT+TT ) = -\ f^gK^KTriSK* " (P,r~ - PK+)) = \l^gK*Kn{£K' ' PK^ 

V 3 V 6 ^^^^ 

A^(p° 7r"7r + ) = J -gpnn{ep ■ {Pn+ ^ Pn')) = " V^^'ptttt (^p " Ptt" ) 

Taking into account all the spin and isospin couplings, the QPCM prediction is 



"'^^For a more detailed definition of the hadronic tensor T^"^ and its parametrization in terms of two 
form factors, fv and /ly, see Ref. [3]. 
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where J(^^^) can be defined from Eq. (71) 

Now, doing a matching between two approaches and factorizing out the common factor 
e-Pi, we can write the following equations: 




6 (76) 

Since J(^^^) is a positive function, one can see that 

s\gn{gK*KTv) = -sign{gp^^) 

and in the SU{3) limit = -^/f- 

One can notice that the choice of the order in the isospin factors of the vector meson 
decay into two pseudoscalars in Eq. (72) well fixes the relative sign of the gvpp couphngs. 
Moreover, this method makes the calculation of the quasi-two-body decay amplitude 
independent on the intermediate vector meson state {K*, p) wave function sign (which, 
in principle, can be arbitrary in the quasi-two-body calculation since the final state is not 
the same)! 



C Partial Wave Amplitudes 

With the quark models one can directly calculate the amplitudes with definite spin or he- 
licity states. An experiment can measure the partial wave amplitudes of particular quan- 
tum numbers of the final state. Since both canonical (orbital) and helicity approaches 
give complete description of the process, one can find the relation between two represen- 
tations for the decay of the initial at-rest state | J, M) with spin J and spin projection M 
on to the z-axis into two particles with spins si^2, helicities Ai^2, total spin 5* and relative 
orbital momentum L [43]: 

Mif,^{Q,) = Njfi,^Dil,^_,,^iQ^) (77) 



with the normalization factor Nj = y 

The observed number of events is given by 



E f Mi'SM)Mi^lrmdnr = An Yl ^y-LJt', (78) 



"'^^One has to be careful with the choice of the momentum, i.e. pc or ps = ~PCi since it changes the 
sign of the P-wave ampHtude. 
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The recoupling from the canonical to the hehcity representation is 

Njfix, = -y^L + liL, 0; S, Ai - A2I J, Ai - A2)(si, Ai; S2, -X2\S, Ai - As)^^ (79) 



L,S 



The two-body decay of the axial-vector meson into vector and pseudoscalar mesons 
can proceed in S and D-waves. Using J = 1, Ai = Ay, A2 = 0, the helicity amplitudes in 
the A reference frame can be written in terms of partial wave amplitudes: 



^i/vo= E v/2LTl(L,0;l,Ay|l,Ay)Ai (80) 

L=0,2 

Setting kv along ^-direction (i.e. 9v = 0), the helicity amplitudes are 

Ml^o = ^1/00 = ^5 - V2Ad 

1+11 1 (81) 

M'^^;, = NJi,^, = As + -^An 

By-turn, the partial wave amplitudes are related to the helicity amplitudes as following: 

As^\{2Mll + M'o^) 

Summing over the final and averaging over the initial spin states, the partial width is then 
given by 

r(^ ^ VP) = {\Asf + \ADf)PS2 (83) 

For the l^-decay into two pseudoscalar mesons Pi and P2 in the P-wave the decay 
amplitude will be given by 

Ml^iQi) = iVi/Jo^j^,o(^i) (84) 

where the helicity amplitude is Ni/qq = ^/Sap. 

Correspondingly, averaging over the l/-spin states, the partial width is then given by 

r(V^ ^ P1P2) = \ApfPS2 (85) 
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D Phase space convention 

The non-relativistic partial width is given by 

r{A ^ BC) = 2n\M'-Z^lc\^ x ^'^f''^ (86) 
where two-body non invariant phase space can be written as 

P^f = / d'kBd'kcS\kB + kc)5(EB + Ec - m^) = An^^^^ (87) 
J ruA 

Since QPCM is in principle a non-relativistic model and we are using the relativistic 

Lorentz-invariant tensor formalism to describe B — > Ki'j decay, one has to make some 

kind of continuation. In order to do that one has to 

• Use relativistic kinematics (i.e. — ki + ml). 

• Use relativistic Breit-Wigner forms. 

• Make the non-relativistic decay amplitudes to be "relativistic" correcting the phase 
space: 

\.[A^ BL) ~ ——^\Ma_^bc\ -87r — \Ma_,bc\ 

from where one immediately obtains the relation between the amplitudes 



— * 

Here Ei, ki are the energies and momentum in the A-reference frame. 
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